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Series Editors’ Foreword 



The series Advances in Industrial Control aims to report and encourage tech- 
nology transfer in control engineering. The rapid development of control tech- 
nology has an impact on all areas of the control discipline. New theory, new 
controllers, actuators, sensors, new industrial processes, computer methods, 
new applications, new philosophies, new challenges. Much of this development 
work resides in industrial reports, feasibility study papers and the reports 
of advanced collaborative projects. The series offers an opportunity for re- 
searchers to present an extended exposition of such new work in all aspects 
of industrial control for wider and rapid dissemination. 

The piezoelectric effect was discovered in some naturally occurring ma- 
terials in the 1880s. However it was not until the Second World War that 
man-made polycrystalline ceramic materials were produced that also showed 
piezoelectric properties. Quartz and other natural crystals found application 
in microphones, accelerometers and ultrasonic transducers, whilst the advent 
of the man-made piezoelectric materials widened the field of applications to 
include sonar, hydrophones, and piezo- ignition systems. Today, piezoelectric 
applications include smart materials for vibration control, aerospace and as- 
tronautical applications of flexible surfaces and structures and novel applica- 
tions for vibration reduction in sports equipment (tennis racquets, skis and 
snowboards), In this entry to the Advances in Industrial Control monograph 
series, Reza Moheimani and Andrew Fleming present a comprehensive tour 
of the present state of the art in the theory and applications of the piezoelec- 
tric property pertaining to vibration control and damping. It is such a well- 
thought-out presentation that we feel sure the volume will become a classic 
for the field. 

The opening chapters of the volume (Chapters 1 and 2) present interesting 
overview material and then the fundamentals of piezoelectricity. It is here 
we learn that piezoelectric transducers can be used as actuators to generate 
moment or as sensors to measure strain. An important application is the 
control of a flexible beam and the model for this is found in Chapter 2. 




x Series Editors’ Foreword 

Fundamental theory on the feedback control aspects of piezoelectric sys- 
tems is found in Chapters 3 to 5; it is useful to note that each new con- 
trol method in these chapters is supported by experimental verifications and 
demonstrations. This makes for very enlightening reading of the sequence of 
increasingly sophisticated control methods used to solve a problem where the 
process models are of high order and the process modes are lightly damped. 

The practical issues of the instrumentation used in piezoelectric systems 
are dealt with in Chapter 6. Here, emphasis is given to the fact that the quality 
of the control obtained is highly dependent on the quality of the instrumen- 
tation used in piezoelectric systems. This chapter reports the authors’ direct 
experience of constructing the instrumentation to implement their control de- 
sign work and is a valuable source of “know-how” for use by other researchers 
in the field. 

A key technique used in piezoelectric vibration control systems is that 
of piezoelectric shunt damping. This is where a piezoelectric transducer, 
which is bounded to the flexible structure, is connected to an electrical shunt 
impedance, Z(s). In Chapter 5, Moheimani and Fleming give this design prob- 
lem a novel twist by reinterpreting the shunt damping structure as a feedback 
structure, opening the way for the use of control design methods in this prob- 
lem. Chapters 7 to 10 exploit the feedback interpretation in a number of ways. 
The spatially distributed nature of flexible structures leads naturally to multi- 
variable problems using multi-port shunts and a generic feedback architecture 
is ideal for this problem (Chapter 7). To reduce the sensitivity of the shunt 
circuit design to changes in piezoelectric capacitance and structure resonance 
frequency, online adaptation is used (Chapter 8). As an alternative to online 
adaptation, negative capacitor shunt impedance designs are investigated in 
Chapter 9. Standard feedback control design tools (such as LQG, and 
H 2 ) are used to compute optimal MIMO shunt damping designs in Chapter 
10 . 

To complete the wide-ranging yet thorough contents of the volume, Chap- 
ter 11 looks at hysteresis effects in piezoelectric transducers and Chapter 12 
is concerned with the important and fairly recent application of piezoelectric 
tube scanners in scanning tunnelling microscopes. As with almost all the other 
chapters in the volume, the discussions in Chapters 11 and 12 are supported 
by experimental findings and results. 

This volume will be an essential reference and source book for all con- 
trol academics, researchers and engineers working in the field of piezoelectric 
control systems. It will be equally attractive to those wishing to gain an in- 
troduction to the area and to those wishing to find examples to illustrate the 
use of control in an advanced applications field. 



M.J. Grimble and M.A. Johnson 
Industrial Control Centre 
Glasgow, Scotland, U.K. 




Preface 



Flexible mechanical systems experience undesirable vibration in response to 
environmental and operational forces. The very existence of vibrations can 
limit the accuracy of sensitive instruments or cause significant error in appli- 
cations where high-precision positioning is essential. In many scientific and en- 
gineering applications control of vibrations is a necessity. Piezoelectric trans- 
ducers have been used in countless applications as sensors, actuators, or both. 
When traditional passive vibration control techniques fail to meet the require- 
ments, piezoelectric transducers, in conjunction with feedback controllers, can 
be used to suppress vibrations in an effective way. 

The ability of piezoelectric materials to transform mechanical energy into 
electrical energy and vice versa has been known in scientific circles for well over 
a century. Although a transducer with such capacity was perceived to be ideal 
for vibration control applications, naturally occurring piezoelectric materials 
were found to be of little use for such purposes. The discovery of piezoceramic 
materials, after the second world war, and invention of methods for their high- 
volume manufacturing enabled researchers to utilize piezoelectric transducers 
in a variety of applications. In particular, over the past decade there has 
been significant interest in using piezoelectric transducers as actuators and 
sensors in applications such as vibration control of flexible structures, micro- 
and nano-positioning systems, medical instrumentation, and micro electro- 
mechanical systems. 

The purpose of this book is to present, in a unified fashion, some recent de- 
velopments in the feedback control of vibration using embedded piezoelectric 
sensors and actuators. The book is intended for researchers, graduate stu- 
dents and practicing engineers with an interest in vibration control systems. 
The book covers various ways in which active vibration control systems can 
be designed for piezoelectric laminated structures, paying special attention to 
how such control systems can be implemented in real time. The text contains 
numerous examples and experimental results obtained from laboratory-scale 
apparatus with details of how similar setups can be built. 




Preface 



This book documents research performed in the Laboratory for Dynam- 
ics and Control of Smart Structures at the School of Electrical Engineering 
and Computer Science, University of Newcastle, Australia. The stimulating 
research environment at the School along with its solid academic culture and 
long tradition of scholarship was an ideal setting for the development of this 
monograph. In the production of this book, the authors wish to acknowledge 
the support they have received from the Australian Research Council (ARC), 
the University of Newcastle and the ARC Centre for Complex Dynamic Sys- 
tems and Control (CDSC). Furthermore, the authors wish to acknowledge 
the contributions made by their colleagues, Sam Behrens, Benjamin Vautier, 
Bharath Bhikkaji, Dominik Niederberger and Manfred Morari for the research 
that underlies parts of the material presented in this book. The authors are 
particularly grateful to Ian Petersen, Hemanshu Pota and Dominik Nieder- 
berger who read the book and made valuable suggestions and comments. 



Newcastle, Australia S. O. R. Moheimani 

December 2005 A. J. Fleming 




Contents 



1 Introduction 1 

1.1 Piezoelectric Transducers 1 

1.2 Vibration Control 2 

1.3 Piezoelectric Shunt Damping 3 

1.4 Hysteresis 5 

1.5 Applications 6 

2 Fundamentals of Piezoelectricity 9 

2.1 Introduction 9 

2.2 History of Piezoelectricity 10 

2.3 Piezoelectric Ceramics 11 

2.4 Piezoelectric Constitutive Equations 13 

2.5 Piezoelectric Coefficients 18 

2.5.1 Piezoelectric Constant dij 18 

2.5.2 Piezoelectric Constant gij 19 

2.5.3 Elastic Compliance Sij 20 

2.5.4 Dielectric Coefficient, 20 

2.5.5 Piezoelectric Coupling Coefficient kij 20 

2.6 Piezoelectric Sensor 22 

2.7 Piezoelectric Actuator 23 

2.8 Piezoelectric 2D Actuation 26 

2.9 Dynamics of a Piezoelectric Laminate Beam 29 

2.10 Active and Macro Fiber Composite Transducers 33 

3 Feedback Control of Structural Vibration 37 

3.1 Introduction 37 

3.2 Structural Properties of Resonant Systems 38 

3.3 Modeling and System Identification 40 

3.3.1 Analytic Modeling 40 

3.3.2 Finite Element Analysis 41 

3.3.3 System Identification 41 




XIV 



Contents 



3.4 Velocity Feedback 43 

3.5 Resonant Controllers 45 

3.6 Positive Position Feedback 48 

3.7 Experimental Implementation of PPF Control on an Active 

Structure 52 

3.7.1 Experimental Setup 52 

3.7.2 System Identification 54 

3.7.3 PPF Controller Design 57 

3.8 Self-sensing Techniques 66 

4 Piezoelectric Shunt Damping 73 

4.1 Introduction 73 

4.2 Passive Shunt Damping 74 

4.2.1 Passivity Defined 74 

4.2.2 Linear Techniques 76 

4.2.3 Resonant Single-mode Shunt Circuits 76 

4.2.4 Resonant Multi-mode Shunt Circuits 78 

4.2.5 Non-linear Techniques 83 

4.3 Active Shunt Damping 84 

4.4 Implementation of Resonant Shunt Circuits 85 

4.4.1 Virtual Circuits 85 

4.4.2 Synthetic Admittance 86 

4.5 Experimental Demonstration 89 

5 Feedback Structure of Piezoelectric Shunt Damping 

Systems 93 

5.1 Introduction 93 

5.2 Feedback Interpretations 94 

5.3 Feedback Structure of Passive Shunts 97 

5.4 Reduction of Inductance Requirements by Adding a Parallel 

Capacitor 100 

5.5 Resonant Shunts 103 

5.6 Properties of Resonant Shunts 106 

5.7 Experimental Implementation of Resonant Shunts 107 

5.8 Hybrid Control 116 

6 Instrumentation 119 

6.1 Introduction 119 

6.2 Strain Voltage Measurement 119 

6.3 Voltage Amplifiers 122 

6.3.1 Linear Amplifiers 122 

6.3.2 Switched- mode Implementation 124 

6.4 Current Amplifiers 125 

6.5 Charge Amplifiers 127 

6.6 Synthetic Admittance 130 




Contents 



xv 



6.7 Switched- mode Synthetic Admittance 131 

6.7.1 Device Operation 132 

6.7.2 Boost Configuration 132 

6.7.3 Efficiency 133 

6.7.4 Practical Advantages and Considerations 134 

6.7.5 Experimental Application 136 

6.8 Signal Processing 137 

6.8.1 Impedance Transformations 139 

6.8.2 Admittance Transformations 140 

6.8.3 Example: Digital Synthesis 142 

6.8.4 Example: Analog Synthesis 142 

7 Multi-port Shunts 145 

7.1 Introduction 145 

7.2 Multi-port Piezoelectric Shunt 146 

7.3 Stability of the Shunted System 149 

7.4 Multivariable Shunts 152 

7.5 Decentralized Shunts 154 

7.6 Experimental Results 155 

7.6.1 Experimental Setup 155 

7.6.2 System Identification 157 

7.6.3 Implementation of a Multi-port Synthetic Admittance . 159 

7.6.4 Implementing the Admittance Transfer Function 160 

7.6.5 Results 160 

8 Adaptive Shunt Damping 165 

8.1 Introduction 165 

8.2 Adaptation Law 166 

8.2.1 Adaptive Single-mode RL Shunt 166 

8.2.2 Adaptive Resonant Multi- mode Shunts 168 

8.3 Experiments 170 

8.3.1 Implementation 170 

8.3.2 Test Structure 170 

8.3.3 Two- mode Damping 172 

8.3.4 Four-mode Damping 173 

9 Negative Capacitor Shunt Impedances 177 

9.1 Introduction 177 

9.2 Negative Capacitor Shunt Controllers 178 

9.3 Implementation 180 

9.3.1 Negative Impedance Converter 180 

9.3.2 Synthetic Admittance 181 

9.4 Experimental Application 181 

9.4.1 Control of a Plate Structure 181 

9.4.2 Control of a Beam Structure 183 




XVI 



Contents 



10 Optimal Shunt Synthesis 187 

10.1 Introduction 187 

10.2 Abstracted Model of Shunted Systems 188 

10.3 Experimental Apparatus 190 

10.4 Instrumentation Electronics 192 

10.5 Active S-impedance Shunt Design 193 

10.5.1 Modeling and Parameter Identification 193 

10.5.2 Hqq S-impedance Design 196 



11 



Dealing with Hysteresis 

11.1 Introduction 

11.2 Hysteresis 

11.3 Charge Control versus Voltage Control 

11.4 Resonant Controllers for Charge-driven Piezoelectric 

Actuators 

11.5 Experimental Implementation of a Multivariable Resonant 

Controller 

11.5.1 The Hysteresis Effect 

11.5.2 State-space Model of the Composite System 

11.5.3 Structure of the State-space Model 

11.5.4 System Identification 

11.5.5 Controller Design 

11.6 Experimental Results 

11.7 Some Observations 



203 

203 

204 
206 

211 

212 

213 

215 

216 
217 
220 
222 
223 



12 Nanopositioning 229 

12.1 Introduction 229 

12.2 Scanning Probe Microscopes 229 

12.3 Piezoelectric Tube Scanners 233 

12.4 Shunt Circuit Modeling 237 

12.4.1 Open-loop 237 

12.4.2 Shunt Damping 239 

12.5 Implementation 245 

12.5.1 Grounded- load Charge Amplifier 245 

12.6 Experimental Results 247 

12.6.1 Tube Dynamics 248 

12.6.2 Amplifier Performance 249 

12.6.3 Shunt Damping Performance 250 



References 



257 



Index 



269 




1 



Introduction 



1.1 Piezoelectric Transducers 

Piezoelectric transducers have become increasingly popular in vibration con- 
trol applications. They are used as sensors and as actuators in structural 
vibration control systems. They provide excellent actuation and sensing capa- 
bilities. The ability of piezoelectric materials to transform mechanical energy 
into electrical energy and vice versa was discovered over a century ago by 
Pierre and Jacque Curie. These French scientists discovered a class of mate- 
rials that when pressured, generate electrical charge, and when placed inside 
an electric field, strain mechanically. 

Piezoelectricity, which literally means “electricity generated from pres- 
sure” is found naturally in many mono crystalline materials, such as quartz, 
tourmaline, topaz and Rochelle salt. However, these materials are generally 
not suitable as actuators for vibration control applications. Instead, man-made 
polycrystalline ceramic materials, such as lead zirconate titanate (PZT), can 
be processed to exhibit significant piezoelectric properties. PZT ceramics are 
relatively easy to produce, and exhibit strong coupling between mechanical 
and electrical domains. This enables them to produce comparatively large 
forces or displacements from relatively small applied voltages, or vice versa. 
Consequently, they are the most widely utilized material in manufacturing of 
piezoelectric transducers. 

Piezoelectric transducers are available in many forms and shapes [28]. The 
most widely used piezoelectric transducers are in the form of thin sheets that 
can be bonded to or embedded in composite structures. As actuators they are 
mainly used to generate moment in flexible structures, while as sensors they 
are used to measure strain. Piezoelectric actuators are also available in the 
form of “stacks” , where many layers of materials and electrodes are assembled 
together. These stacks generate large forces but small displacements in the 
direction normal to the top and bottom surfaces. 

Piezoelectric transducers are used in many applications such as structural 
vibration control [11, 14, 73, 161, 176], precision positioning [124, 75, 183], 
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aerospace systems [4, 106, 43], and more recently they have been critical in 
advancing research in nanotechnology [26, 165, 166]. 

Properties of piezoelectric transducers are studied in Chapter 2 of this 
book. This chapter describes piezoelectric constitutive equations, explains 
piezoelectric coefficients that are of significance, describes the functions of 
piezoelectric sensing and actuation, and finally derives the dynamics of a flex- 
ible structure with piezoelectric sensors and actuators. Specific properties of 
such systems are studied in this chapter. 



1.2 Vibration Control 

Piezoelectric transducers have been extensively used in structural vibration 
control applications. Their wide utilization in this specific application can be 
attributed to their excellent actuation and sensing abilities which stems from 
their high electro-mechanical coupling coefficient, as well as their non-intrusive 
nature. 

For vibration control purposes, piezoelectric transducers are bonded to 
or embedded in a composite structure. The type of structures which lend 
themselves to piezoelectric transducers are generally flexible in nature. The 
transfer functions of these systems are of high order, and their poles are very 
lightly damped. Control problems associated with these systems are by no 
means trivial. 

Control and observation spillover can potentially destabilize a closed-loop 
system if the controller is only designed for a small number of in-bandwidth 
modes and then implemented on the full-order system [12, 13]. A number of 
techniques have been proposed to address this issue, the most promising of 
which is based on using compatible collocated sensors and actuators [98]. 

The velocity feedback controller has well-known stability properties when 
implemented on a flexible structure with collocated actuators and sensors 
[12, 13, 77]. In particular, although this control is susceptible to spill-over, 
existence of spill-over does not destabilize the closed-loop system. The fore- 
most difficulty with velocity feedback control is the high controller bandwidth 
and unnecessary control effort at frequencies far away from the resonance 
frequencies of the base structure. The controller is often rolled-off at higher 
frequencies by appending extra low-pass dynamics. Special care must be ex- 
ercised in adding extra dynamics to the controller since the existence of out- 
of-bandwidth modes of the structure may result in instabilities [77] . 

Resonant controllers approximate a velocity feedback controller at and 
close to the resonance frequencies of a structure [136]. Unlike velocity feed- 
back controllers their control effort at other frequencies is negligible. Further- 
more, they are biproper controllers which enjoy similar stability properties 
to velocity feedback controllers. These controllers are introduced and thor- 
oughly studied in Chapter 3. Relatively similar control techniques have been 
proposed in the literature (see [161] and references therein). However, they 
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are mostly tailored to applications that involve vibration control of a single 
structural mode and do not extend to multivariable systems. The presenta- 
tion in Chapter 3 addresses the problem in its most general form. The same 
control structures are used in the remainder of the book to construct sensor- 
less controllers and to derive controllers for charge-driven actuators, all in a 
multivariable setting. 

Positive position feedback (PPF) controllers, originally proposed in [58] 
have been extensively used in structural vibration control applications [162, 
176, 85, 183, 32, 69, 113, 52, 57]. These controllers have a structure reminiscent 
of the system which is being controlled. Compared to resonant controllers they 
enjoy a quick roll off outside of the controlled bandwidth and they are capable 
of achieving higher levels of performance. 

Stability conditions for PPF controlled systems, as derived in [58], do not 
allow for the effect of out-of-bandwidth modes of the base structure. Conse- 
quently, existence of these modes may destabilize the closed-loop system once 
the controller is implemented on the structure. Having said this, the natural 
roll off of this controller at higher frequencies can often be used to alleviate 
the effect of out-of-bandwidth modes. 

Recently it has been demonstrated that the effect of out-of-bandwidth 
modes of a structure can be captured by adding a feed-through term to the 
truncated model of the system [34, 134, 133, 138]. The corrected model is 
a close representation of the true system in the controlled bandwidth. The 
PPF stability condition, as derived in [58], does not allow for the presence 
of a feed-through term in the structure model. These stability conditions are 
re-derived in Chapter 3 and the results are consistent with those reported in 
[58] when a feed-through term is not present in the structure model. 

All control techniques introduced thus far are devised for structures with 
bonded collocated piezoelectric actuators and sensors. In these systems a 
piezoelectric transducer performs a single task of actuation, or sensing. These 
transducers, however, can function simultaneously as actuators and sensors. 
This property has resulted in emergence of self-sensing circuits which integrate 
the functions of sensing and actuation in the same piezoelectric transducer 
[54, 9]. If designed properly, such circuits can reduce the number of requisite 
piezoelectric transducers by half. However, there are problems associated with 
this technique which are discussed in Section 3.8. 



1.3 Piezoelectric Shunt Damping 

Piezoelectric shunt damping is a popular method for vibration suppression 
in flexible structures. The technique is characterized by the connection of an 
electrical impedance to a structurally bonded piezoelectric transducer. Such 
methods do not require an external sensor, and if designed properly, may 
guarantee stability of the shunted system. 
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Although first appearing in [68] , the concept of piezoelectric shunt damp- 
ing is mainly attributed to Hagood and von Flotow [82], who demonstrated 
that a series inductor-resistor network can significantly reduce vibration of a 
single structural mode. In this technique the network, together with the in- 
herent capacitance of the piezoelectric transducer, is tuned to the resonance 
frequency of the mode which is to be damped. Similar to a tuned mechani- 
cal absorber [95], additional dynamics introduced by the shunt circuit act to 
increase the effective structural damping [82] . 

Single-mode damping can be applied to reduce vibration of several struc- 
tural modes with the use of as many piezoelectric transducers and damping 
circuits. However, in many cases, this may not be a practical solution since a 
large number of transducers will be needed if a large number of modes are to 
be shunt damped. This has encouraged researchers to develop multiple mode 
shunt damping circuits which use only one piezoelectric transducer. 

The first multi-mode shunt circuit is credited to Hollkamp [89], who was 
able to successfully suppress the second and third modes of a cantilever beam 
by 19 and 12 dB respectively. The circuit requires as many parallel branches as 
there are modes to control, and the author proposed a numerical optimization 
to determine suitable component values. The resulting optimization problem is 
highly nonlinear, nonconvex and relatively difficult to solve for a large number 
of modes. 

Another multi- mode circuit, proposed by Wu and co-authors, is made of 
current blocking networks inserted in each parallel branch [198, 200, 199, 203]. 
Inductive and capacitive components of the proposed circuit can be found in 
a straightforward manner. A systematic approach for determining effective 
resistance values is presented in reference [16]. Even in its simplest form [203], 
the complexity and order of current-blocking circuits restrict their use to a 
maximum of three modes. 

An alternative multi-mode shunt circuit was recently introduced in [17]. 
The idea is based on inserting current-flowing LC shunts in each parallel 
branch of the multi- mode shunt, effectively isolating branches from one an- 
other at each resonance frequency of the base structure. Besides greatly sim- 
plifying the tuning procedure current-flowing shunt circuits require less com- 
ponents and are effective in damping a large number of modes. A dual of this 
circuit has also been proposed in [67]. 

Multi-mode shunts discussed thus far are to some extent direct extensions 
of the original single-mode circuits with additional circuitry that enables the 
same shunt to damp several modes of a structure. A new paradigm for the 
design of piezoelectric shunt damping circuits was presented in [140]. By view- 
ing the electrical shunt impedance as parameterizing an equivalent collocated 
strain feedback controller, a shunt impedance can be found working back- 
wards from an effective feedback controller. In particular, starting from reso- 
nant controllers as discussed above, two new passive resonant shunt circuits 
are proposed in [140] and [137]. The proposed shunts are low in order, easy 
to tune, and suitable for modally dense systems. The extension of this tech- 
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nique to multivariable, multi-mode shunt design has been reported in [141]. 
Chapters 5 and 7 cover these topics in a unified way. 

All of the shunts introduced so far are passive in nature. This implies 
that they can be realized using passive components, i.e. resistors, induc- 
tors and capacitors. However, the requisite inductors are generally so large 
that they are often implemented using Gyrator circuits [164]. A more prac- 
tical approach is to simulate the presence of an impedance using a synthetic 
impedance/admittance circuit as described in [59, 60], and articulated in 
Chapters 4 and 6. 

Once effective tools for the implementation of a given impedance are de- 
veloped, it is no longer necessary to restrict the shunt impedance to a passive 
circuit. By allowing the shunt impedance to be active, it is possible to de- 
sign shunt damping circuits that are capable of delivering higher levels of 
performance in terms of adding damping to the system. Furthermore, hav- 
ing identified the feedback structure of the shunted system, it is feasible to 
cast the problem of piezoelectric shunt design as a standard feedback control 
design problem. Therefore, advanced control design techniques such as LQG, 
tfoo, # 2 , etc. can be employed to design high-performance shunts. Chapter 10 
describes a number of ways that the problem of active shunt design can be 
formulated, and discusses the design of an active shunt for a cantilever 
beam. Further experimental results are reported in [66]. 



1.4 Hysteresis 

The assumption made thus far has been that piezoelectric transducers are lin- 
ear devices that transduce energy between mechanical and electrical domains, 
where the coupling mechanism is governed via a linear relationship described 
by the IEEE Standards on Piezoelectricity [154]. While this is generally true at 
lower drives when driven by a voltage amplifier, at higher drives piezoelectric 
actuators display nonlinear behavior known as hysteresis. 

The existence of hysteresis in piezoelectric materials is generally attributed 
to residual misalignment of crystal grains in the poled ceramic [74, 94]. Ex- 
istence of hysteresis has been shown to have an adverse effect on the stabil- 
ity and closed-loop performance of voltage-controlled piezoelectric actuators 
[123]. A large number of techniques have been developed that are aimed at 
reducing the hysteresis associated with volt age- driven piezoelectric actuators, 
e.g. see [71], and references therein. In particular, methods such as inversion- 
based Preisach modeling [127] and phase control [47] are two examples of the 
proposed techniques. 

It has been argued that hysteresis is an electrical property of piezoelectric 
materials, which mainly exists between the applied electric field and the result- 
ing electrical charge [78]. Indeed, it has been demonstrated that by controlling 
electrical charge, or current rather than the applied voltage, the hysteresis ef- 
fect can be substantially reduced [147] . The above argument clearly suggests 
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that when possible, piezoelectric actuators should be driven by charge, or cur- 
rent sources as this significantly reduces the hysteretic nature of the actuator. 

Even though this approach has been known for some time, it has not been 
widely used due to the perceived difficulty of driving highly capacitive loads 
such as piezoelectric actuators. The main problem being the existence of offset 
voltages in the charge or current source circuit, which will eventually charge 
up the capacitive load. This will then distort the control signal being applied 
to the piezoelectric load. This issue has been pointed out by a number of 
authors [122, 124]. Recent research [62] proposes a new structure for charge 
and current sources capable of regulating the DC profile of the actuator. Due 
to this development, it is now possible to use electrical charge as the driving 
control signal for piezoelectric actuators in structural control applications. 

Chapter 11 is mainly concerned with control issues related to piezoelectric 
laminated structures driven by charge amplifiers. Resonant controllers that 
were originally developed for voltage-driven piezoelectric actuators, and were 
introduced in Chapter 5 are redesigned to work with charge-driven piezos. 
It is shown that these modified resonant controllers roll off quickly at higher 
frequencies and maintain closed-loop stability of the system in presence of 
out-of-bandwidth vibration modes. 

One of the issues with resonant controllers for voltage-driven piezos is 
the biproper nature of these controllers. Although stability of the closed- loop 
system under resonant controllers is guaranteed, since these controllers do not 
roll off at higher frequencies instabilities may occur due to introduction of a 
phase lag 1 , or other imperfections. 

Given the specific structure of charge-driven piezoelectric laminated struc- 
tures, it is possible to construct resonant controllers in a way that they achieve 
a roll off of 20 dB, or even 40 dB per decade. This makes them comparable to 
PPF controllers. Furthermore, since the underlying system is free of hysteresis 
nonlinearity excellent performance can be achieved. Chapter 11 is based on 
the results reported in [143] and [188]. 



1.5 Applications 

Piezoelectric vibration control has shown promise in a variety of applications 
ranging from consumer and sporting products to satellite and fighter aircraft 
vibration control systems. 

In the consumer products category, a number of companies such as HEAD 
and K2 have invested in high-performance and novelty items such as compos- 
ite piezoelectric tennis racquets, skis, and snowboards [21]. These products 
typically involve the use of a shunted piezoelectric transducer to decrease 
vibration. Benefits include increased user comfort, better handling and per- 
formance. 



1 This may be partially due to the hysteresis 
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Figure 1.1. A PZT instrumented dual-stage actuator. Picture courtesy of Robert 
Evans, Hutchinson Technology. 



The next generation of hard disk drives may also incorporate piezoelectric 
vibration control systems in a number of ways. The most promising tech- 
nology, which is about to be commercialized, appears to be the utilization 
of dual-stage microactuators in hard disk drives to improve bandwidth and 
tracking performance of the servo system. In these actuation mechanisms, the 
read-write head, which is sitting at the end of a small suspension is moved by 
two piezoelectric actuators that are placed near the root of the suspension, 
see Figure 1.1. Therefore, the piezoelectric actuators serve as a secondary 
actuation mechanism, in addition to the voice coil motor. 

A considerable research effort has been undertaken on the structural con- 
trol of military aircraft. In certain modes of flight, buffet loads on wing and 
stabilizer airfoils can result in high levels of vibration. Such vibration can lead 
to mechanical fatigue and reduction of the flight envelope and lift performance 
of an airfoil [172]. Examples include: FA-18 wing, body, and stabilizer control 
[90], and F-15 panel control [196]. A piezoelectric laminated FA-18 stabilizer 
is shown in Figure 1.2. Piezoelectric transducers have also been incorporated 
into helicopter rotor blades for the suppression of lightly damped lag modes 
in hingeless rotors [106]. 

Other noise control applications include: suppression of acoustic radiation 
from underwater submersibles [205], launch vehicle structural and acoustic 
noise mitigation [153, 51], acoustic transmission reduction panels [105, 171], 
and active antenna structures [73]. 

A primary consideration in the design of space structures is the vibration 
experienced during launch. In future, structures incorporating piezoelectric 
transducers may form the basis of lightweight, high performance mechanical 
components for use in space applications [4]. 

A very successful application for piezoelectric actuators is their use in the 
scanning stage of scanning probe microscopes, which are used extensively in 
nanotechnology research. These “piezoelectric tube scanners”, first reported 
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Figure 1.2. A piezoelectric laminate FA- 18 vertical stabilizer. Picture courtesy of 
Dr. Steven Galea, Defence Science and Technology Organization. 



in [25], can be used in positioning applications which require high precision 
(down to sub-nanometer levels), and high bandwidth. Chapter 12 is concerned 
with problems associated with this specific class of piezoelectric actuators 
and illustrates how the control design techniques, along with the developed 
instrumentation discussed in earlier parts of this book, can be employed to 
significantly improve performance of these nanopositioning systems. 
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Fundamentals of Piezoelectricity 



2.1 Introduction 

This chapter is concerned with piezoelectric materials and their properties. We 
begin the chapter with a brief overview of some historical milestones, such as 
the discovery of the piezoelectric effect, the invention of piezoelectric ceramic 
materials, and commercial and military utilization of the technology. We will 
review important properties of piezoelectric ceramic materials and will then 
proceed to a detailed introduction of the piezoelectric constitutive equations. 

The main assumption made in this chapter is that transducers made from 
piezoelectric materials are linear devices whose properties are governed by a 
set of tensor equations. This is consistent with the IEEE standards of piezo- 
electricity [154]. We will explain the physical meaning of parameters which 
describe the piezoelectric property, and will clarify how these parameters can 
be obtained from a set of simple experiments. 

In this book, piezoelectric transducers are used as sensors and actuators in 
vibration control systems. For this purpose, transducers are bonded to a flexi- 
ble structure and utilized as either a sensors to monitor structural vibrations, 
or as actuators to add damping to the structure. To develop model-based con- 
trollers capable of adding sufficient damping to a structure using piezoelectric 
actuators and sensors it is vital to have models that describe the dynamics of 
such systems with sufficient precision. 

We will explain how the dynamics of a flexible structure with incorporated 
piezoelectric sensors and actuators can be derived starting from physical prin- 
ciples. In particular, we will emphasize the structure of the models that are 
obtained from such an exercise. Knowledge of the model structure is crucial 
to the development of precise models based on measured frequency domain 
data. This will constitute our main approach to obtaining models of systems 
studied throughout this book. 
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2.2 History of Piezoelectricity 

The first scientific publication describing the phenomenon, later termed as 
piezoelectricity, appeared in 1880 [48]. It was co-authored by Pierre and 
Jacques Curie, who were conducting a variety of experiments on a range 
of crystals at the time. In those experiments, they cataloged a number of 
crystals, such as tourmaline, quartz, topaz, cane sugar and Rochelle salt that 
displayed surface charges when they were mechanically stressed. 

In the scientific community of the time, this observation was considered 
as a significant discovery, and the term “piezoelectricity” was coined to ex- 
press this effect. The word “piezo” is a Greek word which means “to press”. 
Therefore, piezoelectricity means electricity generated from pressure - a very 
logical name. This terminology helped distinguish piezoelectricity from the 
other related phenomena of interest at the time; namely, contact electricity 1 
and pyroelectricity 2 . 

The discovery of the direct piezoelectric effect is, therefore, credited to 
the Curie brothers. They did not, however, discover the converse piezoelec- 
tric effect. Rather, it was mathematically predicted from fundamental laws 
of thermodynamics by Lippmann [118] in 1881. Having said this, the Curies 
are recognized for experimental confirmation of the converse effect following 
Lippmann’s work. 

The discovery of piezoelectricity generated significant interest within the 
European scientific community. Subsequently, roughly within 30 years of its 
discovery, and prior to World War I, the study of piezoelectricity was viewed 
as a credible scientific activity. Issues such as reversible exchange of electrical 
and mechanical energy, asymmetric nature of piezoelectric crystals, and the 
use of thermodynamics in describing various aspects of piezoelectricity were 
studied in this period. 

The first serious application for piezoelectric materials appeared during 
World War I. This work is credited to Paul Langevin and his co-workers 
in France, who built an ultrasonic submarine detector. The transducer they 
built was made of a mosaic of thin quartz crystals that was glued between two 
steel plates in a way that the composite system had a resonance frequency 
of 50 KHz. The device was used to transmit a high-frequency chirp signal 
into the water and to measure the depth by timing the return echo. Their 
invention, however, was not perfected until the end of the war. 

Following their successful use in sonar transducers, and between the 
two World Wars, piezoelectric crystals were employed in many applications. 
Quartz crystals were used in the development of frequency stabilizers for 
vacuum-tube oscillators. Ultrasonic transducers manufactured from piezoelec- 
tric crystals were used for measurement of material properties. Many of the 
classic piezoelectric applications that we are familiar with, applications such 



1 Static electricity generated by friction 

2 Electricity generated from crystals, when heated 
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as microphones, accelerometers, ultrasonic transducers, etc., were developed 
and commercialized in this period. 

Development of piezoceramic materials during and after World War II 
helped revolutionize this field. During World War II, significant research was 
performed in the United States and other countries such as Japan and the 
former Soviet Union which was aimed at the development of materials with 
very high dielectric constants for the construction of capacitors. Piezoceramic 
materials were discovered as a result of these activities, and a number of 
methods for their high-volume manufacturing were devised. The ability to 
build new piezoelectric devices by tailoring a material to a specific application 
resulted in a number of developments, and inventions such as: powerful sonars, 
piezo ignition systems, sensitive hydrophones and ceramic phono cartridges, 
to name a few. 



2.3 Piezoelectric Ceramics 

A piezoelectric ceramic is a mass of perovskite crystals. Each crystal is com- 
posed of a small, tetravalent metal ion placed inside a lattice of larger divalent 
metal ions and O 2 , as shown in Figure 2.1. 

To prepare a piezoelectric ceramic, fine powders of the component metal 
oxides are mixed in specific proportions. This mixture is then heated to form 
a uniform powder. The powder is then mixed with an organic binder and is 
formed into specific shapes, e.g. discs, rods, plates, etc. These elements are 
then heated for a specific time, and under a predetermined temperature. As a 
result of this process the powder particles sinter and the material forms a dense 
crystalline structure. The elements are then cooled and, if needed, trimmed 
into specific shapes. Finally, electrodes are applied to the appropriate surfaces 
of the structure. 

Above a critical temperature, known as the “Curie temperature” , each per- 
ovskite crystal in the heated ceramic element exhibits a simple cubic symmetry 
with no dipole moment, as demonstrated in Figure 2.1. However, at tempera- 
tures below the Curie temperature each crystal has tetragonal symmetry and, 
associated with that, a dipole moment. Adjoining dipoles form regions of local 
alignment called “domains” . This alignment gives a net dipole moment to the 
domain, and thus a net polarization. As demonstrated in Figure 2.2 (a), the 
direction of polarization among neighboring domains is random. Subsequently, 
the ceramic element has no overall polarization. 

The domains in a ceramic element are aligned by exposing the element to 
a strong, DC electric field, usually at a temperature slightly below the Curie 
temperature (Figure 2.2 (b)). This is referred to as the “poling process”. 
After the poling treatment, domains most nearly aligned with the electric 
field expand at the expense of domains that are not aligned with the field, 
and the element expands in the direction of the field. When the electric field is 
removed most of the dipoles are locked into a configuration of near alignment 
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Oxygen O Ti,Zr, 



Figure 2.1. Crystalline structure of a piezoelectric ceramic, before and after polar- 
ization 



(Figure 2.2 (c)). The element now has a permanent polarization, the remnant 
polarization, and is permanently elongated. The increase in the length of the 
element, however, is very small, usually within the micrometer range. 

Properties of a poled piezoelectric ceramic element can be explained by 
the series of images in Figure 2.3. Mechanical compression or tension on the 
element changes the dipole moment associated with that element. This cre- 
ates a voltage. Compression along the direction of polarization, or tension 
perpendicular to the direction of polarization, generates voltage of the same 
polarity as the poling voltage (Figure 2.3 (b)). Tension along the direction 
of polarization, or compression perpendicular to that direction, generates a 
voltage with polarity opposite to that of the poling voltage (Figure 2.3 (c)). 
When operating in this mode, the device is being used as a sensor. That is, 
the ceramic element converts the mechanical energy of compression or tension 
into electrical energy. Values for compressive stress and the voltage (or field 





Figure 2.2. Poling process: (a) Prior to polarization polar domains are oriented 
randomly; (b) A very large DC electric field is used for polarization; (c) After the 
DC field is removed, the remnant polarization remains. 
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Figure 2.3. Reaction of a poled piezoelectric element to applied stimuli 



strength) generated by applying stress to a piezoelectric ceramic element are 
linearly proportional, up to a specific stress, which depends on the material 
properties. The same is true for applied voltage and generated strain 3 . 

If a voltage of the same polarity as the poling voltage is applied to a ceramic 
element, in the direction of the poling voltage, the element will lengthen and 
its diameter will become smaller (Figure 2.3 (d)). If a voltage of polarity 
opposite to that of the poling voltage is applied, the element will become 
shorter and broader (Figure 2.3 (e)). If an alternating voltage is applied to 
the device, the element will expand and contract cyclically, at the frequency 
of the applied voltage. When operated in this mode, the piezoelectric ceramic 
is used as an actuator. That is, electrical energy is converted into mechanical 
energy. 



2.4 Piezoelectric Constitutive Equations 

In this section we introduce the equations which describe electromechanical 
properties of piezoelectric materials. The presentation is based on the IEEE 
standard for piezoelectricity [154] which is widely accepted as being a good 
representation of piezoelectric material properties. The IEEE standard as- 
sumes that piezoelectric materials are linear. It turns out that at low electric 
fields and at low mechanical stress levels piezoelectric materials have a linear 
profile. However, they may show considerable nonlinearity if operated under 
a high electric field or high mechanical stress level. In this book we are mainly 
concerned with the linear behavior of piezoelectric materials. That is, for the 
most part, we assume that the piezoelectric transducers are being operated 
at low electric field levels and under low mechanical stress. 

When a poled piezoelectric ceramic is mechanically strained it becomes 
electrically polarized, producing an electric charge on the surface of the mate- 
rial. This property is referred to as the “direct piezoelectric effect” and is the 

3 It should be stressed that this statement is true when the piezoelectric material is 
being operated under small electric field, or mechanical stress. When subject to 
higher mechanical, or electrical fields, piezoelectric transducers display hysteresis- 
type nonlinearity. For the most part, in this monograph, the linear behavior of 
piezoelectric transducers will be of interest. However, Chapter 11 will briefly re- 
view the issues arising when a piezoelectric transducer is operated in the nonlinear 
regime. 
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Dipole Alignment 



Figure 2.4. Schematic diagram of a piezoelectric transducer 



basis upon which the piezoelectric materials are used as sensors. Furthermore, 
if electrodes are attached to the surfaces of the material, the generated elec- 
tric charge can be collected and used. This property is particularly utilized in 
piezoelectric shunt damping applications to be discussed in Chapter 4. 

The constitutive equations describing the piezoelectric property are based 
on the assumption that the total strain in the transducer is the sum of me- 
chanical strain induced by the mechanical stress and the controllable actuation 
strain caused by the applied electric voltage. The axes are identified by numer- 
als rather than letters. In Figure 2.4, 1 refers to the x axis, 2 corresponds to 
the y axis, and 3 corresponds to the 2 axis. Axis 3 is assigned to the direction 
of the initial polarization of the piezoceramic, and axes 1 and 2 lie in the plane 
perpendicular to axis 3. This is demonstrated more clearly in Figure 2.5. 

The describing electromechanical equations for a linear piezoelectric ma- 
terial can be written as [154, 70]: 

£i = Sfjtjj + d mi E m (2.1) 

D m = + ^ k Ek, (2-2) 

where the indexes = 1,2,. ..,6 and m,k = 1,2,3 refer to different di- 
rections within the material coordinate system, as shown in Figure 2.5. The 
above equations can be re-written in the following form, which is often used 
for applications that involve sensing: 



— ^ij^j ~b 9 mi Dm 

Ei = 9mi®i T fiikDp, 



where 



(2.3) 

(2.4) 
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# Axis 



1 x 

2 y 

3 £ 

4 Shear around x 

5 Shear around y 

6 Shear around £ 



Figure 2.5. Axis nomenclature 



a .. . stress vector ( N/m 2 ) 

6 . . . strain vector {m/m) 

E . . . vector of applied electric field (V/m) 

£ . . . permitivity ( F/m ) 

d ... matrix of piezoelectric strain constants (m/V) 

S . . . matrix of compliance coefficients (m 2 /N) 

D . . . vector of electric displacement (C/m 2 ) 
g . . . matrix of piezoelectric constants (m 2 /C) 

/3 . . . impermitivity component (m/F) 

Furthermore, the superscripts D , E , and a represent measurements taken 
at constant electric displacement, constant electric field and constant stress. 

Equations (2.1) and (2.3) express the converse piezoelectric effect, which 
describe the situation when the device is being used as an actuator. Equations 
(2.2) and (2.4), on the other hand, express the direct piezoelectric effect, which 
deals with the case when the transducer is being used as a sensor. The converse 
effect is often used to determine the piezoelectric coefficients. 

In matrix form, Equations (2.1)- (2. 4) can be written as: 







'Su S 12 S13 514 Sib S 16 ' 




<71 


£2 




S 2 1 S22 S23 S2A S 25 $26 




°2 


£3 




5 3 l S32 S33 S34 S35 S36 




(73 


£4 




S41 ^42 S43 S44 S45 S4e 




723 


£5 




Sei S 52 £53 ^54 S 55 See 




731 


£ 6 _ 




_Se 1 Se2 Se 3 Se 4 See See_ 




Jl2 _ 



+ 



du d<ii dsi 

dl2 d22 ds2 

dl3 d23 d33 
dl4 d24 ^34 
<^15 ^25 ^35 
die ^26 ^36 



Ei 

E 2 

E3 



(2.5) 



and 
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'd r 




dn di2 dis di4 di5 diQ 


D 2 


= 


^21 ^22 ^23 ^24 ^25 ^26 


D 3 




_^31 ds2 ds3 ds4 ^35 d^Q_ 





1 

b2 

03 

03 

b^n 

<13 




'e r 


+ 


e?l e 22 ^23 




^2 




p<j p<j 

L e 31 e 32 e 33j 




^3 



Some texts use the following notation for shear strain 



723 = £4 

731 = £5 
712 = £6 



and for shear stress 



( 2 . 6 ) 



723 = 04 
731 = CT5 
712 = CT6- 



Assuming that the device is poled along the axis 3, and viewing the piezo- 
electric material as a transversely isotropic material, which is true for piezo- 
electric ceramics, many of the parameters in the above matrices will be either 
zero, or can be expressed in terms of other parameters. In particular, the 
non-zero compliance coefficients are: 



Su = S 22 

Sis = S31 = S 23 = S 32 
S12 = S21 

S44 = S55 

S 66 = 2(5 n -Si 2 ). 

The non-zero piezoelectric strain constants are 

^31 = ^32 

and 

^15 = ^24- 

Finally, the non-zero dielectric coefficients are and e% 3 . Subsequently, 

the equations (2.5) and (2.6) are simplified to: 




and 
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£1 




£2 




£3 




e 4 
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Sn S12 Sis 00 0 
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0 0 0 0 644 0 
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dis 
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dis 
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0 
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°2 




^3 




723 




T 3 I 




_Tl2_ 



(2.7) 



~D 1 - 
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0 
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D 2 
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dis 
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D 3 




_dsi 


dsi 


ds3 
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0 



e n 


0 


1 

0 




0 


pCT 

e ll 


0 




_ 0 


0 


pC T 

e 33 j 





Ei 

E 2 

Es 



0 % 

CT2 
cr 3 
cr 4 

^5 

^6 



The “piezoelectric strain constant” d is defined as the ratio of developed 
free strain to the applied electric field. The subscript dij implies that the elec- 
tric field is applied or charge is collected in the i direction for a displacement 
or force in the j direction. The physical meaning of these, as well as other 
piezoelectric constants, will be explained in the following section. 

The actuation matrix in (2.5) applies to PZT materials. For actuators 
made of PVDF materials, this matrix should be modified to 



0 0 dsi 

0 0 ds2 

0 0 dss 

0 <^25 0 

di5 0 0 

0 0 0 



This reflects the fact that in PVDF films the induced strain is nonisotropic 
on the surface of the film. Hence, an electric field applied in the direction of 
the polarization vector will result in different strains in 1 and 2 directions. 
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Figure 2.6. A piezoelectric transducer arrangement for cfai measurement 



2.5 Piezoelectric Coefficients 

This section reviews the physical meaning of some of the piezoelectric coeffi- 
cients introduced in the previous section. Namely dij, gij, Sij and e^-. 



2.5.1 Piezoelectric Constant dij 



The piezoelectric coefficient dij is the ratio of the strain in the j-axis to the 
electric field applied along the i-axis, when all external stresses are held con- 
stant. In Figure 2.6, a voltage of V is applied to a piezoelectric transducer 
which is polarized in direction 3. This voltage generates the electric field 




which strains the transducer. In particular 




in which 



At 



d 31 V£ 

t 



The piezoelectric constant d 3 i is usually a negative number. This is due 
to the fact that application of a positive electric field will generate a positive 
strain in direction 3. 

Another interpretation of dij is the ratio of short circuit charge per unit 
area flowing between connected electrodes perpendicular to the j direction to 
the stress applied in the i direction. As shown in Figure 2.7, once a force F is 
applied to the transducer, in the 3 direction, it generates the stress 



F 



= 



iw 



which results in the electric charge 



q = d 33 F 
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F 




Figure 2.7. Charge deposition on a piezoelectric transducer - An equal, but opposite 
force, F, is not shown 



flowing through the short circuit. 

If a stress is applied equally in 1, 2 and 3 directions, and the electrodes 
are perpendicular to axis 3, the resulting short-circuit charge (per unit area), 
divided by the applied stressed is denoted by d p . 



2.5.2 Piezoelectric Constant gij 



The piezoelectric constant gij signifies the electric field developed along the 
i-axis when the material is stressed along the j-axis. Therefore, in Figure 2.8 
the applied force F, results in the voltage 



w 



Another interpretation of gij is the ratio of strain developed along the 
j-axis to the charge (per unit area) deposited on electrodes perpendicular to 
the i-axis. Therefore, in Figure 2.9, if an electric charge of Q is deposited on 
the surface electrodes, the thickness of the piezoelectric element will change 
by 



A£ = 



93lQ 

w 




Figure 2.8. An open-circuited piezoelectric transducer under a force in direction 1 
- An equal, but opposite force, F, is not shown 
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Figure 2.9. A piezoelectric transducer subject to applied charge 



2.5.3 Elastic Compliance Sij 

The elastic compliance constant Sij is the ratio of the strain the in ^-direction 
to the stress in the j-direction, given that there is no change of stress along 
the other two directions. Direct strains and stresses are denoted by indices 1 
to 3. Shear strains and stresses are denoted by indices 4 to 6 . Subsequently, 
Si 2 signifies the direct strain in the 1 -axis when the device is stressed along 
the 2 -axis, and stresses along directions 1 and 3 are unchanged. Similarly, S 44 
refers to the shear strain around the 2 -axis due to the shear stress around the 
same axis. 

A superscript “E” is used to state that the elastic compliance Sfj is mea- 
sured with the electrodes short-circuited. Similarly, the superscript “D” in 
SP denotes that the measurements were taken when the electrodes were left 
open-circuited. A mechanical stress results in an electrical response that can 
increase the resultant strain. Therefore, it is natural to expect to be smaller 
than Sfj . That is, a short-circuited piezo has a smaller Young’s modulus of 
elasticity than when it is open-circuited. 

2.5.4 Dielectric Coefficient, 

The dielectric coefficient e^- determines the charge per unit area in the i - axis 
due to an electric field applied in the j-axis. In most piezoelectric materials, 
a field applied along the j-axis causes electric displacement only in that di- 
rection. The relative dielectric constant, defined as the ratio of the absolute 
permitivity of the material by permitivity of free space, is denoted by K. 
The superscript a in e\ x refers to the permitivity for a field applied in the 1 
direction, when the material is not restrained. 

2.5.5 Piezoelectric Coupling Coefficient kij 

The piezoelectric coefficient represents the ability of a piezoceramic mate- 
rial to transform electrical energy to mechanical energy and vice versa. This 
transformation of energy between mechanical and electrical domains is em- 
ployed in both sensors and actuators made from piezoelectric materials. The 
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ij index indicates that the stress, or strain is in the direction j, and the 
electrodes are perpendicular to the i-axis. For example, if a piezoceramic is 
mechanically strained in direction 1, as a result of electrical energy input in 
direction 3, while the device is under no external stress, then the ratio of 
stored mechanical energy to the applied electrical energy is denoted as k^ x . 

There are a number of ways that kij can be measured. One possibility 
is to apply a force to the piezoelectric element, while leaving its terminals 
open-circuited. The piezoelectric device will deflect, similar to a spring. This 
deflection A z , can be measured and the mechanical work done by the applied 
force F can be determined 

FA Z 

Due to the piezoelectric effect, electric charges will be accumulated on the 
transducer’s electrodes. This amounts to the electrical energy 



wy = 



20 , 



which is stored in the piezoelectric capacitor. Therefore, 



&33 = 



Wm 

Q 



y/F^C~ p 



The coupling coefficient can be written in terms of other piezoelectric 
constants. In particular 



( 2 . 8 ) 

where E p is the Young’s modulus of elasticity of the piezoelectric material. 

When a force is applied to a piezoelectric transducer, depending on 
whether the device is open-circuited or short-circuited, one should expect to 
observe different stiffnesses. In particular, if the electrodes are short-circuited, 
the device will appear to be “less stiff” . This is due to the fact that upon the 
application of a force, the electric charges of opposite polarities accumulated 
on the electrodes cancel each other. Subsequently no electrical energy can 
be stored in the piezoelectric capacitor. Denoting short-circuit stiffness and 
open-circuit stiffness respectively as K sc and K oc , it can be proved that 
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QE p cr 

U ij ij 

— 9ijdijE p , 
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2.6 Piezoelectric Sensor 



When a piezoelectric transducer is mechanically stressed, it generates a volt- 
age. This phenomenon is governed by the direct piezoelectric effect ( 2 . 2 ). 
This property makes piezoelectric transducers suitable for sensing applica- 
tions. Compared to strain gauges, piezoelectric sensors offer superior signal 
to noise ratio, and better high-frequency noise rejection. Piezoelectric sen- 
sors are, therefore, quite suitable for applications that involve measuring low 
strain levels. They are compact, easy to embed and require moderate signal 
conditioning circuitry. 

If a PZT sensor is subject to a stress field, assuming the applied electric 
field is zero, the resulting electrical displacement vector is: 
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The generated charge can be determined from 



Q = 



j J [Di D 2 D 3 ] 



dAi 
dA 2 , 
dA s 



where dAi, dA 2 and dAs are, respectively, the differential electrode areas in 
the 2 - 3 , 1-3 and 1-2 planes. The generated voltage V p is related to the charge 



where C p is capacitance of the piezoelectric sensor. 

Having measured the voltage, V p , strain can be determined by solving the 
above integral. If the sensor is a PZT patch with two faces coated with thin 
electrode layers, e.g. the patch in Figure 2 . 4 , and if the stress field only exists 
along the 1-axis, the capacitance can be determined from 



C p m 



l we 33 

t 



Assuming the resulting strain is along the 1-axis, the sensor voltage is found 
to be 

d 3 iE p w f 

V s = — -7 / £i dx, ( 2 . 9 ) 

^ p Ji 

where E p is the Young’s modulus of the sensor and s\ is averaged over the 
sensor’s length. The strain can then be calculated from 
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CpVs 

d 3 \E p iw 



( 2 . 10 ) 



In deriving the above equation, the main assumption was that the sensor was 
strained only along 1-axis. If this assumption is violated, which is often the 
case, then (2.10) should be modified to 

C p Vs 

£l (1 — v)d 3 \E p iw ’ 



where v is the Poisson’s ratio 4 . 



2.7 Piezoelectric Actuator 



Consider a beam with a pair of collocated piezoelectric transducers bonded to 
it as shown in Figure 2.10. The purpose of actuators is to generate bending 
in the beam by applying a moment to it. This is done by applying equal 
voltages, of 180° phase difference, to the two patches. Therefore, when one 
patch expands, the other contracts. Due to the phase difference between the 
voltages applied to the two actuators, only pure bending of the beam will 
occur, without any excitation of longitudinal waves. The analysis presented 
in this section follows the research reported in references [42, 11, 76, 53]. 

When a voltage V is applied to one of the piezoelectric elements, in the 
direction of the polarization vector, the actuator strains in direction 1 (the 
x-axis). Furthermore, the amount of free strain is given by 



_ d 31 V 

£ p — . * 

ip 



( 2 . 11 ) 



where t p represents the thickness of the piezoelectric actuator. 

Since the piezoelectric patch is bonded to the beam, its movements are 
constrained by the stiffness of the beam. In the foregoing analysis perfect 
bonding of the actuator to the beam is assumed. In other words, the shearing 
effect of the non-ideal bonding layer is ignored [33] . Assuming that the strain 
distribution is linear across the thickness of the beam 5 , we may write 



e(z) = az. (2.12) 

The above equation represents the strain distribution throughout the 
beam, and the piezoelectric patches, if the composite structure were bent, 
say by an external load, into a downward curvature. Subsequently, the por- 
tion of the beam above the neutral axis and the top patch would be placed in 

4 Notice that if cfei / e.g. if the sensor is a PVDF film, then this expression 

for strain must be changed to £i = 5 . 

(1 - v $UL)d 31 E p lw 

5 This is consistent with the Kirchoff hypothesis of laminate plate theory [97]. 
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Figure 2.10. A beam with a pair of identical collocated piezoelectric actuators 



tension, and the bottom half of the structure and the bottom patch in com- 
pression. Although, the strain is continuous on the beam-actuator surface, 
the stress distribution is discontinuous. In particular, using Hooke’s law, the 
stress distribution within the beam is found to be 



a b {z) = E h az , (2.13) 

where E b is the Young’s modulus of elasticity of the beam. Since the two 
“identical” piezoelectric actuators are constrained by the beam, stress distri- 
butions inside the top and the bottom actuators can be written in terms of 
the total strain in each actuator (the strain that produces stress) 

(7 l v = E p (az - e p ) (2-14) 

(Tp = E p (az + £p ), (2.15) 

where E p is the Young’s modulus of elasticity of the piezoelectric material 
and the superscripts t and b refer to the top and bottom piezoelectric patches 
respectively. Applying moment equilibrium about the center of the beam 6 
results in 




(Tp(z)zdz + 




c jp{z)zdz + 




c Tp{z)zdz = 0. 



After integration a is determined to be 



(2.16) 



3^((|+^) 2 -(|) 2 ) 

a 2(E p {Ci+t P r-ar}+E b Cir) £p - 

6 Due to the symmetrical nature of the stress field, the integration need only 
be carried out starting from the centre of the beam, i.e. f Q 2 a p (z)zdz + 

/tf +tp crp(z)zdz = 0. 

2 
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Figure 2.11. A beam with a single piezoelectric actuator 



The induced moment intensity 7 , M in the beam is then determined by 
integrating the triangular stress distribution across the beam: 

M = E b Ia , (2.18) 

where / is the beam’s moment of inertia. Knowledge of M is crucial in deter- 
mining the dynamics of the piezoelectric laminate beam. 

If only one piezoelectric actuator is bonded to the beam, such as shown in 
Figure 2.11, then the strain distribution (2.12) needs to be modified to 

e(z) = (az + eo)- (2.19) 

This expression for strain distribution across the beam thickness can be 
decomposed into two parts: the flexural component, az and the longitudinal 
component, £q. Therefore, the beam extends and bends at the same time. This 
is demonstrated in Figure 2.12. The stress distribution inside the piezoelectric 
actuator is found to be 



a p (z) = E p (az + £ 0 - £ p ). (2.20) 

The two parameters, So and a can be determined by applying the moment 
equilibrium about the centre of the beam 




Figure 2.12. Decomposition of asymmetric stress distribution (a) into two parts: 
(b) flexural and (c) longitudinal components. 

moment per unit length 



7 
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/ a b (z)zdz + / cr p (z)zdz = 0 



and the force equilibrium along the x-axis 



( 2 . 21 ) 



r^-+tp 

cr b (z)dz + / <j p {z)dz = 0. (2.22) 

Unlike the symmetric case, the force equilibrium condition (2.22) needs to 
be applied. This is due to the asymmetric distribution of strain throughout 
the beam. Solving (2.21) and (2.22) for so and a we obtain 




6E b E p t b tp(t b + t p ) 

E 2 h t\ + E p E b (4tp p + 6tgi* + 4^3) + E 2 tp £ P 



(2.23) 



and 

{E b tp + Eptp}E p (t b / 2 ) 

£ ° = Wt + E p E b (41^t p + 6 ipf + 4t 6 t3) + E%t p £p ‘ 



(2.24) 



The response of the beam to this form of actuation consists of a moment 
distribution 

M x = E h Ia (2.25) 

and a longitudinal strain distribution 



£x — £o- 



(2.26) 



It can be observed that the moment exerted on the beam by one actuator 
is not exactly half of that applied by two collocated piezoelectric actuators 
driven by 180° out-of-phase voltages. This arises from the fact that the Ex- 
pressions (2.25) and (2.23) do not include the effect of the second piezoelectric 
actuator. However, if this effect is included by allowing for the stiffness of the 
second actuator in the derivations, while ensuring that the voltage applied to 
this patch is set to zero, then it can be shown that the resulting moment will 
be exactly half of that predicted by (2.18) and (2.17). The collocated situa- 
tion is often used in vibration control applications, in which one piezoelectric 
transducer is used as an actuator while the other one is used as a sensor. This 
configuration is appealing for feedback control applications for reasons that 
will be explained in Chapter 3. 



2.8 Piezoelectric 2D Actuation 

This section is concerned with the use of piezoelectric actuators for excitation 
of two-dimensional structures, such as plates in pure bending. The analysis 
is similar to that presented in the previous section. A typical application is 
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Figure 2.13. A piezoelectric actuator bonded to a plate 



shown in Figure 2.13, which demonstrates a piezoelectric transducer bonded 
to the surface of a plate. It is also assumed that another identical transducer 
is bonded to the opposite side of the structure in a collocated fashion. If the 
two patches are driven by signals that are 180° out of phase, the resulting 
strain distribution, across the plate, will be linear as shown in Figure 2.14 a 
and b. That is, 



= a x z (2.27) 

Ey = OtyZ , (2.28) 

where a x and a y represent the strain distribution slopes in the x — z and y — z 
planes respectively. 

Assuming that the piezoelectric material has similar properties in the 1 
and 2 directions, i.e. dsi = ^ 32 , the unconstrained strain associated with the 
actuator in both the x and y directions, under the voltage V, is given by 

d 3 iV 



Now the resulting stresses in the plate, in the x and y directions are 

E , 

Ox = x _ v 2 y £x + ve v) 

and 

E , 

°y = Y~^^ SyJrVSx ^ 

where v is the Poisson’s ratio of the plate material. Representing the stresses 
in the top piezoelectric patch as erg and cr^, and the stresses in the bottom 
patch as <rg and cr^, we may write 
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(a) 




Figure 2.14. Two dimensional strain distribution in a plane with two collocated 
anti-symmetric piezoelectric actuators 



a x — -t _ 2 l/ p) £ p} (2.29) 

i Vp 

°x — . _ P 2 "L 6 * + V P £ y + (■*■ + v p) £ p) (2.30) 

1 l/p 

a y = i — {^2/ + — (1 + (2.31) 

I — 

E 

& y = ^ 2 {^2/ ^p^x T (1 T z/p)^p} , (2.32) 



where z/ p is the Poisson’s ratio of the piezoelectric material. 

Given that e p is the same in both directions and that the plate is homo- 
geneous, we may write 

£ x — £y — S. 

Subsequently, the strain distribution across the plate thickness can be written 
as 

£ = OL x Z = OLyZ = az. 
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The condition of moment of equilibrium about the x and y axes can now be 
applied. That is, 

t t | 4 - 

r 2 r 2 

J < 7 x zdz + J a^zdz = 0 



and 



J (jyzdz + J ^ 






a v y zdz = 0 , 



where £ represents the plate’s thickness. Integrating and solving for a gives 
3£ p {(| +^) 2 -(|) 2 }(l-Q 



= 



£r>. 



+ i P ) 3 - (^) 3 }(1 - 0 + 2 £(f ) 3 (1 - v P T P 



The resulting moments in x and y directions are 



M x = My = Ela. 



(2.33) 



For the symmetric case, i.e. when only one piezoelectric actuator is bonded 
to the plate, similar derivations to the previous section can be made. 



2.9 Dynamics of a Piezoelectric Laminate Beam 



In this section we explain how the dynamics of a beam with a number of 
collocated piezoelectric actuator /sensor pairs can be derived. At this stage we 
do not make any specific assumptions about the boundary conditions since we 
wish to keep the discussion as general as possible. However, we will explain 
how the effect of boundary conditions can be incorporated into the model. 

Let us consider a setup as shown in Figure 2.15, where m identical collo- 
cated piezoelectric actuator /sensor pairs are bonded to a beam. The assump- 
tion that all piezoelectric transducers are identical is only adopted to simplify 
the derivations, and can be removed if necessary. The actuator is exposed 
to a voltage of v ai (t) and the voltage induced in the i ^ sensor is v Pi (t). 

We assume that the beam has a length of L, width of W , and thickness 
of £ 5 . Corresponding dimensions of each piezoelectric transducer are LF P , 
and t p . Furthermore, we denote the transverse deflection of the beam at point 
x and time £ by z(x,t). The dynamics of such a structure are governed by the 
Bernoulli-Euler partial differential equation 



E b I 



d A z(x , £) 
dx A 



+ pAb 



d 2 z(x , £) 
dt 2 



d 2 M x {x,t) 
dx 2 



(2.34) 



where p, A&, Eb and / represent density, cross-sectional area, Young’s mod- 
ulus of elasticity and moment of inertia about the neutral axis of the beam 
respectively. The total moment acting on the beam is represented by M x (x , £), 
which is the sum of moments exerted on the beam by each actuator, i.e. 
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Figure 2.15. A beam with a number of collocated piezoelectric actuator/sensor 
pairs 



m 

M x (x,t) = y ~^M Xi (x,t). (2.35) 

i= 1 

The moment exerted on the beam by the actuator, M Xi (x,t) can be 
written as 



M Xi (x, t ) = Kv ai (t){u(x — xi.) — u(x — X2 i )} 1 (2.36) 

where u(x) represents the unit step function, i.e. u(x) = 0 for x < 0 and 
u(x) = 1 for x > 0. The term {u(x — x±.) — u(x — X2J} is incorporated into 
(2.36) to account for the spatial placement of the actuator. The constant R 

can be determined from (2.11), (2.17) and (2.18). The forcing term in (2.34) 
can now be determined from Expressions (2.36) and (2.35), and using the 
following property of Dirac delta function 




S^ n \t — 6)(j)(t)dt 



(-!)><"> ( 0 ), 



(2.37) 



where 5^ is the n ^ derivative of 5, and 0 is a continuous function of 0 
[111]. Having determined the expression for the forcing function in (2.34) we 
can now proceed to solving the partial differential equation. One approach to 
solving this PDE is based on using the modal analysis approach [130]. In this 
technique the solution of the PDE is assumed to be of the form 



z(x,t) = y ^W k (x)q k (t). (2.38) 

k = 1 

Here Wk(x ), known as the modeshape, is the eigenfunction which is deter- 
mined from the eigenvalue problem obtained by substituting (2.38) into (2.34) 
and using the following orthogonality properties [130] 




/ wk(x)w p (x)dx = S kp 

J 0 



E b I d A w k (x) 
pA b dx 4 



w p (x)dx 




(2.39) 

(2.40) 
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th 

where ujk describes the k in natural frequency of the beam and Sk p is the 
Kroneker delta function, i.e. 

_ ( i, k =p 
kp \ 0 , otherwise 

A solution to the eigenvalue problem requires precise knowledge of bound- 
ary conditions. For specific boundary conditions, e.g. simply-supported and 
cantilevered, mode shapes and resonance frequencies can be determined ana- 
lytically from the eigenvalue problem. For further details, the reader is referred 
to [130] and [50]. 

A set of uncoupled ordinary differential equations can be obtained from 
(2.34) using the orthogonality properties (2.39) and (2.40) and the property 
(2.37), as well as (2.38). It can be shown that the ordinary differential equa- 
tions are of the form: 

rn 

'52i’k i Va i (t),. (2.41) 

2=1 

where k = 1 , 2 , .. . and qk(t) is the generalized coordinate of the k ^ mode. 
Furthermore, the parameter is found to be: 

i/jki = / Wk(x) {8\x — xi.) — S'(x — X 2 i)} dx (2.42) 

J o 

= w 'k( x 2i) -w f k (xii) y (2.43) 

where f'(x) represents the first derivative of the function / with respect to x, 
and we have used (2.37) to obtain (2.43) from (2.42). 

To this end we point out that the differential equation (2.41) does not 
contain a term to account for the natural damping associated with the beam. 
The presence of damping can be incorporated into (2.41) by adding the term 
2 CkQk(t) to (2.41). This results in the differential equation 

m 

y ^kjVai ( t )■ (2.44) 

2=1 

Applying the Laplace transform to (2.44), assuming zero initial conditions, 
we obtain the following transfer function from the vector of applied actuator 
voltages V a (s) = [v ai (s), . . . , v am (s)] to the beam deflection z(x,s) at loca- 
tion x 



q k {t) + 2 Ck<^kQk(t) + ulq k (t) = 



Qk(t) + u> k q k (t) = 



pAb 



OO 

G(x,s) = 7 ^ 



Wk{xW k 

S 2 + 2 Ck^k s + 



where 7 = -A- and 

1 pA b 



(2.45) 
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i’k = 



> fe i' 






til 

The piezoelectric voltage induced in the i lLL sensor can be obtained using 
Expression (2.9). That is, 



v Pi {t) = 



^3 1 EpWp 



a 



s r ,dx. 



p J o 



t h 

The expression for the mechanical strain in the i lli sensor patch can be 
obtained from 







tb 



+ tr 



d 2 z 
dx 2 ’ 



Now v Pi is found to be 






—dziEpWp (^ + t p ) 



OO 

'Yh'Pk i qk{t). 

k = 1 



(2.46) 



Therefore, the transfer function matrix relating the voltages applied to the 
piezoelectric actuators V a (s) = [^(s), • • . , ^a m (s)] to the voltages measured 
at the piezoelectric sensors V p (s) = [v Pl (s), . . . , v Pm (s)] is found to be: 



OO 

Gw(s) = lE 

k= 1 



Mi 

s 2 + 2 CfeWfeS + a; 2 



(2.47) 



where 

—dsiE p W p (y + tp) R 
C p pA}) 

To this end we need to explain how the systems represented by Infinite 
Series (2.47) and (2.45) can be approximated with finite dimensional models. 
In any controller design scenario we may only be interested in designing a 
controller for a finite bandwidth. If N modes of the structure lie within that 
bandwidth of interest, the series (2.47) and (2.45) are often truncated to obtain 
a finite dimensional model of the structure, which is of minimum dimensions. 
While the truncation may not be of concern in non-collocated models, e.g. 
(2.45), it can be a serious problem for a collocated system, e.g. (2.47). Trun- 
cation of a collocated model can result in perturbations in open-loop zeros of 
the system, which in the worst case can cause closed-loop instabilities, and in 
the best case will contribute to the loss of closed- loop performance [34] . 

A number of techniques have been proposed to compensate for the ef- 
fect of truncated out-of-bandwidth modes on collocated structural models. 
The reader is referred to [34, 207, 138, 135, 84] and references therein for an 
overview of such techniques. We do, however, point out that the truncation 
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Figure 2 . 16 . (a) A flexible structure with a pair of collocated piezoelectric trans- 
ducers, and (b) its electrical equivalent model 



error can be minimized by appending the truncated model by a feed-through 
term, i.e. to approximate (2.47) with 



N 

G vv (s) =7 

k= 1 



S 2 + 2 Ck^kS + 



+ D. 



If chosen properly, the addition of this feed-through term to the truncated 
model can result in an acceptable approximation. In this book, in most cases, 
we will use system identification to directly identify this parameter, as well as 
the rest of the dynamics of the system. 

We conclude this chapter with an important observation that will be uti- 
lized in the forthcoming chapters. Figure 2.16 (a) illustrates a flexible struc- 
ture with a collocated piezoelectric actuator /sensor pair. The piezoelectric 
transducer on the right functions as an actuator, and the one on the left as 
a sensor. What is of importance here is the electrical model of the system 
depicted in Figure 2.16 (b). Each piezoelectric transducer is modeled as a ca- 
pacitor in series with a strain-dependent voltage source. The transfer function 
from the voltage applied to the actuator, v to the voltage induced in the sen- 
sor, v s = v p is given by (2.47). This observation is particularly important in 
designing piezoelectric shunts, and to identify the connection between piezo- 
electric shunt damping and feedback control of a collocated system. This will 
be discussed in more detail in Chapter 5. 



2.10 Active and Macro Fiber Composite Transducers 

Active Fiber Composites (AFCs) are an alternative to traditional monolithic 
piezoelectric transducers. First proposed in 1992 [20], longitudinally polarized 
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Figure 2.17. A piezoelectric Active Fiber Composite (AFC) comprised of piezo- 
electric fibers with interdigitated electrodes top and bottom 



piezoelectric fibers, as shown in Figure 2.17, are encased in an epoxy resin 
with interdigitated electrodes laminated onto the top and bottom surfaces of 
the transducer. With an applied voltage, the interdigitated electrodes induce 
longitudinal electric fields along the length of each fiber. The original motiva- 
tion was to increase the electromechanical coupling by utilizing the high d 33 
piezoelectric strain constant rather than the lesser dsi constant. 

Active fiber composites have a number of practical advantages over tradi- 
tional monolithic transducers [19]: 

• The fibers are encapsulated by the printed polymer electrodes and epoxy 
resin thus increasing the reliability and service-life in harsh environments. 

• The short length and diameter or the fibers together with their alignment 
along the length of the transducer increases the conformability of AFC 
transducers. They can be laminated onto structures with complex geome- 
tries and curvatures. 

• AFC transducers are more robust to mechanical failure than monolithic 
transducers. In addition to their conformability, they can also tolerate local 
and incremental damage. If some of the fibers are fractured, the transducer 
will not be substantially damaged, in contrast, monolithic actuators will 
fracture and fail if they are stressed beyond their yield limit. 

• AFCs have been reported to develop greater strains than monolithic actu- 
ators [19]. The strain actuation is also unidirectional. 

Macro Fiber Composites (MFCs) [174] are similar in nature to AFCs as 
they utilize the direct dss piezoelectric effect through the use of interdigitated 
electrodes. Rather than individual fibers, a monolithic transducer is simply 
cut into a number of long strips. The resulting transducer is conformable in 
one dimension and more robust to mechanical failure than monolithic patches. 

The greatest disadvantages of AFC and MFC transducers is their high 
present cost, and the large voltages required to achieve the same actuation 
strain as monolithic transducers. The equivalent piezoelectric capacitance is 
also much lesser making them unsuitable as low-frequency strain sensors (see 
Section 6 . 2 ). 
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The low capacitance of AFC and MFC transducers also causes difficulties 
in the implementation of piezoelectric shunt damping systems, to be discussed 
in Chapter 4. Device capacitances of less than 50 nF have been deemed im- 
practical for shunt damping [18]. A performance comparison of monolithic, 
AFC, and MFC transducers in a passive shunt damping application can be 
found in references [18] and [148]. 

Although the piezoelectric transducers used throughout this book are ex- 
clusively monolithic, all of the techniques discussed in the proceeding chapters 
are equally as applicable to AFC and MFC variants. Indeed, from the control 
engineers viewpoint, transducer physics is usually lumped into a simplified 
electrical model, or identified as part of the structural system. 
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Feedback Control of Structural Vibration 



3.1 Introduction 

This chapter is concerned with the problem of feedback control system design 
for vibration reduction in flexible structures. In the previous chapter we stud- 
ied some important properties of flexible structures. We paid special attention 
to those structures in which piezoelectric transducers were used in the form 
of actuators and sensors. 

Transfer functions of flexible structures have interesting properties. For 
example, different transfer functions associated with the same structure have 
identical poles. These systems have a large number of poles, and these poles 
are very close to the juj axis. The phase of a collocated system is between 0 
and -180 degrees. This property implies that poles and zeros of the system 
interlace. 

Some characteristics of flexible structures make the task of controlling 
them rather difficult. For example, the fact that these transfer functions are 
of very high orders means that if a controller is designed with a view to control 
only a small number of in-bandwidth modes of the system, the existence of 
out-of-bandwidth modes could potentially destabilize the closed-loop system. 
This problem has been investigated in the literature in the context of spill-over 
effect [12, 13]. 

Despite the fact that control of flexible structures is a difficult task, it 
is possible to design and implement high-performance controllers capable of 
reducing structural vibrations and ensuring closed-loop stability of the con- 
trolled system. In particular, when actuators and sensors are collocated and 
compatible, e.g. collocated piezoelectric transducers, very effective control de- 
sign techniques can be proposed. 

All of the control design techniques introduced in this chapter assume a 
collocated structure. Collocation is a property of flexible structures that en- 
ables the control designer to devise feedback structures that are guaranteed 
to be closed- loop stable in the presence of out-of-bandwidth dynamics. Al- 
though in some applications it may be difficult to collocate the sensor and 
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Figure 3.1. A flexible structure with several collocated piezoelectric actua- 
tor/sensor paris 



actuator, if piezoelectric transducers are to be used for both actuation and 
sensing collocation should not be viewed as a difficulty. Collocated piezoelec- 
tric actuator /sensor packs are commercially available in the market. 

In this chapter we will also introduce the problem of self-sensing piezoelec- 
tric actuators. The idea is to use a piezoelectric transducer simultaneously as a 
sensor and an actuator, thereby reducing the number of requisite piezoelectric 
transducers by half. This is done by simulating the presence of a collocated 
piezoelectric sensor using properties of the piezoelectric actuated system and 
an electronic circuit. Some of the techniques and procedures are described and 
complications associated with the proposed methods are discussed. A more 
practical approach will be introduced in the next chapter. 



3.2 Structural Properties of Resonant Systems 



As explained in the previous chapter resonant systems, such as flexible struc- 
tures are continuous distributed systems described by hyperbolic partial dif- 
ferential equations. However, for practical considerations, it is common to 
approximate these systems by lumped models. Such models can be obtained 
through a number of ways, e.g. finite element modeling, modal analysis or 
system identification. 

In this book we are mainly concerned with flexible structures with embed- 
ded, or bonded piezoelectric transducers. Such a system is shown in Figure 3.1. 
In the previous chapter we demonstrated that the multivariable collocated 
transfer function of the system in Figure 3.1 can be represented as 



M 

g vv (s ) = 

i= 1 



S 2 + 2(iiL>iS + id? 



(3.1) 



where ^ is an m x 1 vector (assuming m piezoelectric actuator /sensor pairs), 
and M — > oo. In practice, however, the integer M is finite, but possibly a very 
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Figure 3.2. Feedback control of vibration with a collocated piezoelectric actua- 
tor/sensor pair 



large number which represents the number of modes that sufficiently describe 
the elastic properties of the structure under excitation [93, 100]. 

Also shown in Figure 3.1 is a number of disturbances acting on the struc- 
ture. These may represent point forces, distributed forces (such as a wind 
gust) or torques acting on the structure. Nevertheless, the transfer function 
matrix relating the disturbance vector W to the vector of voltages measured 
at the sensors V p can be written as 



M 



G vw (s ) — 



Ipil'i 



s 2 + 2 (iU>iS + U)f ’ 



(3.2) 



where ji is a £ x 1 matrix, assuming there are t disturbances acting on the 
structure. 

For the SISO case, i.e. when G vv (s) represents the transfer function asso- 
ciated with only one collocated actuator /sensor pair, this system is known to 
possess interesting properties [126]. In particular, it is known that the system 
is minimum phase, and furthermore, poles and zeros of the system interlace. 
This ensures that phase of the collocated transfer function will be within the 
0 to —180° range. 

Given the highly resonant nature of flexible structures it is natural to 
investigate ways of adding damping to the system. One approach to adding 
damping to the structure is to use a feedback controller, as illustrated in 
Figure 3.2. 

The transfer function matrix (3.1) consists of a very large number of par- 
allel second order terms. This is demonstrated in block diagram form in Fig- 
ure 3.3, which illustrates the situation when such a system is being controlled. 
In most scenarios, only control of a limited bandwidth is of importance. Typ- 
ically N modes of the structure would fit within this bandwidth while modes 
N + 1 and above are left uncontrolled. The uncontrolled modes, however, do 
exist and have the potential to destabilize the closed-loop system. Therefore, 
the existence of these modes should be taken into account, and a controller 
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Controlled Modes 




Figure 3.3. Feedback control of flexible structure consisting of a large number of 
elastic modes 



should be designed to ensure adequate damping performance, as well as sta- 
bility in the presence of these out-of-bandwidth modes. 

A number of control design techniques which satisfy the above require- 
ments will be introduced in the forthcoming sections. 



3.3 Modeling and System Identification 

During the analysis and design of structural vibration control systems it is 
of great benefit to possess a dynamic system model. That is, a state- space or 
transfer-function model describing outputs such as sensor voltages and dis- 
placements in response to external forces and control inputs. A model allows 
the designer to utilize model-based design tools and to simulate the charac- 
teristics of a closed-loop system before experimental implementation. A brief 
discussion follows on the three most common techniques for procuring models 
of flexible structures. It is concluded that the last technique, system identifi- 
cation, is the most suitable for the applications described in this book. 

3.3.1 Analytic Modeling 

Analytic modeling, typically involving the assumed modes approach [130], 
requires distinct models for both structural dynamics and piezoelectric trans- 
ducers [70]. The modal analysis procedure has been used extensively through- 
out the literature for obtaining models of structural [130] and acoustic systems 
[87]. Its major disadvantage is the requirement for detailed physical informa- 
tion regarding the sensors, actuators, and the underlying mechanical system. 
Practical application typically involves the use of experimental data and a 
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non-linear optimization to identify unknown parameters such as modal am- 
plitudes, resonance frequencies and damping ratios. Even in this case, the 
descriptive partial differential equations must still be solved (as functions of 
the unknown parameters) to obtain the mode shapes. This may be difficult 
or impossible for realistic structural or acoustic systems with complicated 
boundary conditions. 

As discussed in Section 2.9, the Lagrangian or modal expansion represents 
a small deflection of a distributed parameter system as an infinite summation 
of modes [130]. The modes are a product of two functions, one of the spatial 
co-ordinate x, and another of the temporal £, i.e. 

oo 

z(x,t) = '%2w k (x)q k (t), (3.3) 

k = 1 

where the qk(t)s are the modal displacements, z(x,t) is the displacement at 
a point, and the Wk(x ) s are the system eigenfunctions or mode-shapes in the 
context of mechanical systems analysis. The mode shapes Wk(x) must form 
a complete coordinate basis for the system, satisfy the geometric boundary 
conditions, and for analytic analysis be differentiable over the spatial domain 
to at least the degree required by the describing partial differential equations. 
Many practical systems also obey certain orthogonality conditions. The reader 
is referred to Section 2.9 for an in-depth discussion of the Partial Differential 
Equations (PDEs) governing flexible structures and their solution by modal 
analysis. 

3.3.2 Finite Element Analysis 

Another popular technique for obtaining structural models is that of finite 
element (FE) analysis [40]. This is an approximate method that results in 
high-order spatially discrete systems. If sensors and actuator dynamics are 
known, the integrated model can be cast in a state-space form to facilitate 
control design and analysis [115]. 

In situations involving non-trivial boundary conditions or irregular geom- 
etry, finite element modeling is generally regarded as a better alternative to 
analytic modeling. Due to the approximate nature of finite element modeling, 
there is no requirement to solve a partial differential equation, thus, the pro- 
cess is significantly more straightforward. Unfortunately, detailed information 
regarding the structures’ material properties and boundary conditions is still 
required. Similar to the modal analysis procedure, errors in these parameters 
are accounted for by correcting the FE models with experimental data [56] . 

3.3.3 System Identification 

System identification can be employed to procure a composite structural and 
piezoelectric model directly from experimental data. Although the field of 
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system identification is extremely diverse [120], the range of possibilities is 
significantly reduced if we restrict ourselves to techniques capable of identify- 
ing systems with multiple inputs and outputs, that do not require an explicit 
parameterization or non-linear optimization, two rather undesirable traits. 
The majority of the residue comprises the so-called subspace class of sys- 
tem identification algorithms. Such methods identify state-space models by 
exploiting geometric properties of the input and output sequences. Methods 
exist for both time and frequency domain data, see [189] for a summary of 
time domain methods. 

Before selecting a system identification algorithm it is worth consider- 
ing the various options for data collection and processing. The traditional 
approach has been to record input and output sequences then apply a time- 
domain system identification technique directly. When considering systems 
with highly resonant modes a number of problems arise: 

1. Unless the system dynamics are known a-priori , it is impossible to excite 
all of the modes of interest equally, or to achieve a reasonably constant 
signal to noise ratio over the entire bandwidth. Typically the first mode 
will dominate the energy of the measurement. 

2. Data records can be prohibitively long. As an example consider a typical 
mechanical system - a cantilever beam as described in Section 8.3.2 with 
a first-mode resonance frequency of 11 Hz and a fourth- mode resonance 
frequency of 400 Hz. Allowing for the roll-off of anti-aliasing and recon- 
struction filters, a suitable sampling frequency would be around 4-10 kHz. 
To record only 10 cycles of the lowest frequency mode, approximately 
10,000 samples would be required. Records of this size typically result 
in numerical instabilities, out-of- memory problems, and large amounts of 
computing time. This problem is exacerbated if the system under consid- 
eration has multiple inputs or outputs. 

Frequency domain data is a more suitable format for highly resonant systems. 
Frequency domain data can be acquired directly by swept sine analysis or by 
computing the ratio of input and output Fourier transforms. A full discus- 
sion on obtaining empirical transfer function estimates (frequency responses) 
from time domain data can be found in [159] and [120]. The most straight- 
forward approach is to excite the system with a periodic disturbance, e.g. 
a periodic chirp, and to record exactly an integer number of periods. This 
avoids Fourier transform leakage and negates the need for window functions 
introducing systematic error. Even though a frequency response obtained by 
Fourier transform still requires relatively long time histories, these can be 
heavily truncated before performing the identification. Furthermore, the data 
points of interest can be chosen arbitrarily, although this is impossible in the 
time domain. 

The best technique for measuring frequency response data from highly 
resonant mechanical systems is by swept sine analysis. A spectrum analyzer 
such as the Hewlett-Packard 35670 A will provide extremely high-quality data 
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with logarithmically spaced frequencies and auto-leveled inputs to achieve an 
excellent signal-to- noise ratio over the entire bandwidth. Each frequency point 
is estimated from a large amount of time domain data, hence the superior noise 
performance. 

Of the various frequency domain system identification algorithms treated 
in [159], subspace methods stand out, they are simple and have proven to be ef- 
fective in identifying resonant systems of high-order [129, 128]. The multivari- 
able discrete-time algorithm of McKelvey [128] is used extensively throughout 
this book. Continuous time systems can be identified by applying a straight- 
forward frequency pre-warping to the data and a bilinear transform to the 
resulting model [128]. Commercial implementations can be found in the Sys- 
tem Identification Toolbox for MATLAB® [119] and the Frequency Domain 
System Identification Toolbox for MATLAB® [158]. Another algorithm used 
in this book is that of Van Overschee and De Moor [186]. This is a continuous 
time method that avoids the usual numerical conditioning problems associated 
with a basis constructed from the powers of s. 

A further benefit of procuring models in the frequency domain is that 
the model quality of wide-band systems can be easily assessed by overlaying 
the data and model response. If a suitable model cannot be obtained it is 
also possible to perform some design and analysis procedures with only fre- 
quency response data. Non-parametric controllers such as resonant shunts can 
be designed with only frequency response data. It is also straightforward to 
estimate the closed-loop response with only frequency response data and a 
controller transfer function. 



3.4 Velocity Feedback 

A controller that guarantees unconditional closed-loop stability of a collocated 
system is the velocity feedback controller. This feedback controller is defined 
as [12, 13, 10] 



AT(s) = q'qs, (3.4) 

where 7 is a 1 x m matrix with m the number of collocated actuator sensor 
pairs. 

Closed- loop stability of (3.1) under (3.4) can be verified in a number of 
ways. However, in the SISO case, this can be explained graphically. Being 
a collocated transfer function, the phase of G vv (s ) is between 0 and -180 
degrees. This, along with the fixed 90 degree phase of the differentiator and 
the negative feedback gain of - 1 , ensures that Nyquist plot of sG vv (s ) will not 
circle -1+jO point in the complex plane. 

To prove closed-loop stability in the multivariable case the following the- 
orem is needed. 
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Theorem 3.1. Consider the following second order multivariable dynamical 
system 

x(t) + Vx(t) + JCx(t) = 0, (3.5) 

where V, JC £ R 7VxAr and x £ H Nxl . Furthermore, assume that V = V > 0. 
Then (3.5) is exponentially stable if and only if JC = JC' >0. 

Proof: The second order system (3.5) can be re-written as a first order system 
by defining X(t) = [x(t) x(t)]: 



X(t) =AX(t), 



where 



Now, let 



where 



A = 



P = 



0 / 

-JC -V 

JC + TV tI 
tI I 



0 < r < A min (P). 

Then, it is straightforward to verify that 



A'P + PA = 

-2 rlC 0 
0 2 (r/ - D) 



< 0 . 



(3.6) 

(3.7) 



(3.8) 



This is a necessary and sufficient condition for exponential stability of 
(3.6). ■ 

The next theorem establishes that (3.4) is a stabilizing controller for (3.1). 

Theorem 3.2. Negative feedback connection of (3.1) and (3.4) is exponen- 
tially stable. 



The following theorem proves closed-loop stability of the multivariable 
system. 

Proof. The transfer function materix G vv (s ) can be represented by the fol- 
lowing second-order state-space equation 



where 



x(t) + 2 ZQx(t) + C2 2 x(t) = F'u(t) 
y(t) = &x(t), 



(3.9) 
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rci 



z = 



(2 



Cm I 



Q = 



CJi 



UJ 2 



% 

^ = [V>1 V>2 • • • V’m] • 

Also, the feedback controller (3.4) incorporating the negative sign can be 
represented as follows 



u(t) = - 7 ; 7 y(t). (3.10) 

Hence, the closed-loop system dynamics can be described by 

x(t) + (2Zi? + iF'y'ylF) x(t) + fi 2 x(t ) = 0. (3.11) 

Clearly, 2Zi? + iF'y'ylF is a symmetric, positive-definite matrix. There- 
fore, according to Theorem 3.5 we can conclude that the closed-loop system 
represented by (3.11) is exponentially stable. 

Although velocity feedback allows for spill-over, existence of spill-over does 
not destabilize the closed-loop system. This is a very interesting property. 
However, despite its appealing stability properties, velocity feedback as de- 
scribed by (3.4) is of little use in the context of piezoelectric laminated struc- 
tures. The main problem is due to the requirement for a differentiator, which is 
not possible. Some extra dynamics would have to be added to the compensator 
to ensure that it will roll off at high frequencies. Due to the high-bandwidth 
nature of the underlying system these additional dynamics have the potential 
to destabilize the closed-loop system. 

Another problem with velocity feedback, which is perhaps of a greater con- 
cern, is high control effort at all frequencies. For vibration damping purposes, 
the control effort should ideally be restricted to regions around the resonance 
frequencies of the base structure. This calls for different classes of controller 
that will be introduced in the remainder of this chapter. 



3.5 Resonant Controllers 

As indicated earlier, the structure of a collocated system, as expressed by 
(3.1), allows for the design of feedback controllers with specific structures that 
guarantee unconditional stability of the closed-loop system. Such controllers 
are of interest due to their ability to avoid closed- loop instabilities arising from 
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the spill-over effect [13]. Resonant controllers have the distinctive property 
that closed-loop stability of the system is guaranteed in presence of out-of- 
bandwidth modes. 

Two possible resonant controllers are: 



and 



N 1 2 
a _ OLi^S 

v “ s 2 + 2diU)iS + uf ’ 
1=1 1 


(3.12) 


js(3 _ + 2 diCbi) 

V 2—1 s 2 + 2diU>iS + L of ’ 

1=1 1 


(3.13) 



where typically N <C M. Note that both cq and are m x 1 vectors. Fur- 
thermore, Co\ < < . • . < Cbfj. 

The typical feedback control problem associated with the system (3.1) 
and a resonant controller K v is illustrated in Figure 3.2. Here, v p is the vector 
of voltages measured at the piezoelectric sensors, while w is the vector of 
disturbances acting on it. The purpose of the controller is to add damping 
to the structure, hence reducing the effect of disturbances. This is done by 
shifting closed-loop poles of the system deeper into the left half of the complex 
plane. 

Closed- loop stability of the multivariable collocated system (3.1) under 
(3.12) and (3.13) can be proved in a number of ways. We will first prove 
stability of the feedback connection of (3.1) and (3.13). 



Theorem 3.3. Feedback connection of (3.1) and (3.13) is stable. 
Proof: Stability of (3.1) and (3.13) is equivalent to the stability of 



M 

g vv (s) = 

i= 1 



s 

s 2 + 2Q(JiS + caf 



and 

~/3 _ + 2djQi) 

v “ S 2 + 2diCbiS + 

i=i 1 

Now, it can be proved that G vv (s ) is a positive-real (PR) transfer function 
matrix, while K/( is a strictly positive-real (SPR) transfer function matrix. 
Negative feedback connection of a PR and a SPR system is known to be 
stable [104]. ■ 

Closed- loop stability of (3.1) under (3.12) can be established by construct- 
ing a Lyapunov function. 

Theorem 3.4. Feedback connection of (3.1) and (3.12) is exponentially sta- 
ble. 
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Figure 3.4. Equivalent system for study of closed- loop stability of K% and G vv 



Proof: Feedback connection of (3.1) and (3.12) is exponentially stable if and 
only if the closed-loop system depicted in Figure 3.4 is stable, where 



M 



Gw(s) = 






■Mis 

S + “IQiLOiS + w? 



and 



K« = 



N , 

OLiOL^S 

r-f s 2 + 2 diihiS + u>f ' 
1=1 1 



Now, G vv (s) can be represented by 



x(t) + 2 Zf2x(t) + fi 2 x(t) = 

y(t) = 



where 




Q 



Cn_ 

Ml 



CJN 



& = [Vh ^2 • • • H • 



Furthermore, K£(s) can be written as 1 

x(t) + 2 AQx(t) + fi 2 x(t) = r'y(t ) 
?x(£) = —Tx{t ), 



(3.14) 

(3.15) 



(3.16) 

(3.17) 



where 



1 Note that the negative sign is incorporated in the feedback controller. 
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u 



A = 



d/vl 



Q = 



CJi 



CJ2 



Un\ 

r = [ai c^2 • • • • 

Combining (3. 14)- (3. 15), the closed- loop system dynamics can be obtained 

(3.18) 



'x{t) 


1 


'2Zfi wr 


~x(t) 


1 


'ti 1 0 ' 


x{t) 




'o' 




“T 


-r'v An_ 


_x(t)_ 


“T 


. 0 ^ 2 . 






0 



Now, a Lyapunov function, V(x(t),x(t)) can be defined as 
V(x(t),x(t)) = 



'±(ty 


/ 




i 


'x{t) 


/ 


so 

to 


0 1 


'x(t) 


i(t)_ 




M*). 


T 


x{t)_ 




0 


to 


_x(t)_ 



(3.19) 



Clearly, V(x(t),x(t)) > 0 for all nontrivial x, x, x and x. Differentiating V 
with respect to time yields 



V(x(t),x(t)) 



'x{t) 


r 


'-2 ZQ '//'/' 


'x(ty 


x{t)_ 




r'W -2AQ 


_x(t)_ 



(3.20) 



which implies 

V(x(t),x(t)) < 0 



for nontrivial values of x and x. Therefore the closed- loop system is asymp- 
totically stable. ■ 



3.6 Positive Position Feedback 

One of the difficulties associated with the implementation of resonant con- 
trollers is that the frequency responses of (3.12) and (3.13) do not roll off at 
higher frequencies. Although under ideal assumptions the closed-loop stability 
of the collocated system (3.1) under (3.12) and (3.13) is guaranteed, in real 
world applications due to issues such as phase contributions of the antialias- 
ing filters, etc. the existence of out-of-bandwidth dynamics may destabilize 
the closed-loop system. It is, therefore, of interest to have a controller which 
possesses a similar structure to (3.12) and (3.13), but whose frequency re- 
sponse rolls off at higher frequencies. Positive position feedback, proposed by 
Caughey and co-authors [77, 58] has such properties. 
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For a system of the form (3.1), a positive position feedback controller is 
defined as 



Kpp(s) 



N 

S' 



-lili 



■ 25i£jiS + lj- 



2 ’ 



(3.21) 



where 7 i G R mxl and Si > 0 for i— 1, 2, . . . , N. 

Due to the existence of the negative sign in all terms of (3.21) the overall 
system resembles a positive feedback loop. Also, the transfer function matrix 
(3.1) is similar to that of the force to displacement transfer function matrix 
associated with a flexible structure, hence the terminology positive position 
feedback. 

To derive stability conditions for this control loop, the series in (3.1) is 
first truncated by keeping the first N modes ( N < M ) that lie within the 
bandwidth of interest, and then incorporating the effect of truncated modes 
by adding a feed-through term to the truncated model [34, 135]. That is to 
approximate (3.1) by 



N 



Gvv ( s ) — y 






S 2 + 2 CiiUiS + 



D. 



(3.22) 



The following theorem gives the necessary and sufficient conditions for 
closed-loop stability under positive position feedback. 



Theorem 3.5. The negative feedback connection of (3.22) and (3.21), with 
Z > 0 is stable if and only if 

n 2 - r'Dr > 0 ( 3 . 23 ) 

and 

n 2 - ^'r(fi 2 - r'Dr)~ l r > 0 , ( 3 . 24 ) 

where 



ro 



z = 



(2 



Cam 



Q = 



CJi 



W>2 



CJN m 

& = |>1 fh... ^n] 



and 
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r*i 



A = 



^ N A 



kJi 



Q = 



CjO 2 



uk 



NA 



r = [71 72 • • • 7 n ] • 

Proof: Using the above notation, (3.22) and (3.21) can be represented in 
second-order state-space form as 



and 



x(t) + 2 ZQx(t) + f2 2 x(t) = \P'u(t) 

y(t ) = ^x(t) + Du(t ) 

x(t) + 2 AQx(t) + fi 2 x(t) = r'y(t) 



u(t) = rx{t). 

The closed-loop system can then be written as 






2 zn 0 


'x(t) 




^0 

to 

1 


'x(t) 




'o' 


0 2 An_ 




+ 


-r'r n 2 - r'Dr\ 


x(t)_ 




0 



(3.25) 

(3.26) 



(3.27) 

(3.28) 



.(3.29) 



Clearly 



2Zi? 0 
0 2Z\i2 



>0, 



therefore, according to Theorem 3.1 the closed- loop system is stable, if and 
only if 

Q 2 -VT 

-r'r Q 2 - r'Dr 



> 0. 



(3.30) 



The Schur complement [27] implies that this inequality holds if and only if 

n 2 - r'Dr > o 



and 

q 2 - r’r(n 2 - r'Dry x r’r > o 

which proves the theorem. ■ 

It should be pointed out that the stability condition (3.23), (3.24) can 
be represented by a linear matrix inequality in variables A, Q 2 and IF, by 
applying the Schur complement as indicated in the following lemma: 
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Figure 3.5. Beam layout as used in the experiments 



Lemma 3.6. The set of positive position feedback controllers is a convex set 
characterized by the following linear matrix inequalities: 



and 



A> 0 



(3.31) 



fl 2 -&T 0 

-r'v n 2 r 
o r d - 1 



>o 



(3.32) 



Proof: This can be proved using the Schur complement as follows: 



ft 2 -XT 

-r'v n 2 - r'Dr 



$ 



' n 2 


-XT' 




0" 


-r'v 


£>< 

to 

1 




r 



d [o r\ > o 



cj 2 -TT 0 

-t't q 2 r 
o r d - 1 



> 0 . 



An implication of Lemma 3.6 is that all admissible positive position feed- 
back controllers must be stable, as asserted by Condition (3.31). 

Positive position feedback controllers are generally capable of ensuring 
higher levels of performance compared with resonant controllers (3.12) and 
(3.13). However, in making a comparison one has to be careful about the 
assumptions made in proving stability of the closed-loop systems associated 
with two different classes of controllers. To prove stability under Feedback 
Controllers (3.12) and (3.13) no simplifying assumption on the system had to 
made. Indeed it was proved that stability of the full order system was guaran- 
teed under these two controllers. However, to obtain stability conditions for 
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Figure 3.6. Picture of the cantilever beam 



the positive position feedback controller, we had to assume that a truncated 
model of the system with only N modes adequately represented the dynamics 
of the system within the bandwidth of interest. Although, from a mathemati- 
cal point of view this may be problematic, in most realistic applications where 
sensors and actuators have limited bandwidths, this assumption is unlikely to 
be of concern. Indeed many applications for positive position feedback con- 
trol have been reported in the literature. To this end the reader is referred to 
[8, 3, 170, 109, 177] and references therein. 

To conclude this section we point out that the stability condition (3.23), 
(3.24) does allow for a non-zero feed through term in the model of the sys- 
tem (3.22). As discussed in Section 2.9, inclusion of a feed through term in 
the model of a flexible structure may be essential. Particularly, when sensors 
and actuators are collocated. It should be noted that this analysis is missing 
from the original work of Caughy and Fanson [58]. Nevertheless, the stability 
condition reported in [58] can be derived from (3.23), (3.24) by setting D to 
zero. 

3.7 Experimental Implementation of PPF Control on an 
Active Structure 

3.7.1 Experimental Setup 

All discussions on resonant systems so far were based on Analytic Models 
(3.1), (3.2) and (3.22). The model structures of the resonant controllers and 
PPF controllers, see (3.12), (3.13) and (3.21), strongly resemble the system 
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model structures (3.1), (3.2) and (3.22). Here, a cantilever beam representing 
a physical resonant system is considered. This beam, which is clamped at one 
end and free at the other end, is susceptible to high amplitude vibrations 
when disturbed. In what follows, a PPF controller will be designed to damp 
these highly resonant vibration modes of the beam. It is worth noting that 
cantilever beams of the type used here are known to have models of the form 
(3.1), (3.2) and (3.22), as articulated in Section 4, and resonant controllers 
designed for them are known to have damped the vibrations to a reasonable 
extent, see [187] and [144]. 

In this section a description of the experimental setup meant for both 
modeling and control of the beam is presented. This setup will remain the 
same for all the experiments to be performed on the beam. 

As mentioned above, the cantilever beam is clamped at one end and free 
at the other. Two pairs of piezoelectric patches are attached to this beam, one 
pair located close to the clamped end and the other pair is located close to 
the free end of the beam. For each pair, one piezoelectric patch will be used 
as an actuator (where input signals are applied) and the other patch will act 
as a sensor (where output signals are recorded) . Another solitary piezoelectric 
patch is attached to the center of the beam and will be driven by a voltage 
source w. This voltage w represents the disturbance acting on the beam. See 
Figures 3.5 and 3.6 for a schematic beam setup and an actual picture of the 
beam respectively. 

Not surprisingly, the mechanical properties (including the damping and 
the resonance frequencies) will depend on the dimensions and material prop- 
erties of the beam and the piezoelectric patches. In Figure 3.7, Table 3.1 and 
Table 3.2 the dimensions of the beam setup and the material properties of the 
beam and the piezoelectric transducers are given. 

NOT TO SCALE 




The thickness of each collocated patch is 0.25 mm and the “disturbance” patch is 0.5 mm 



Figure 3.7. Beam setup and dimensions 
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Table 3.1. Beam Properties 



Length, L 


550 mm 


Thickness, h 


3 mm 


Width, W 


50 mm 


Density, p 


2.77 x 10 3 kg/m 3 


Young’s Mod., E 


7.00 x 10 10 N/m 2 



Table 3.2. PIC 151 Ceramic Properties 



Length, L pz 


50 mm 


Thickness, h pz 


0.25 mm 


Width, W pz 


25 mm 


Charge Constant, cfei 


-210 x 10“ 12 m/V 


Voltage Constant, g 31 


-11.5 x 10" 3 Vm/N 


Coupling Coefficient, &31 


0.34 


Capacitance, C p 


115 nF 



3.7.2 System Identification 

In order to obtain a model for the beam, the experimental setup is treated as 
a three-input-three-output multivariable system, see Figure 3.8. The inputs 
(v\ and V 2 ) in Figure 3.8 are the voltages applied to the actuators of the 
collocated piezo patches and the outputs v p \ and v P 2 are the voltages induced 
at the corresponding sensors. The third input w is the disturbance on the 
beam and the output y t i p is the displacement of the the tip of the beam. To 
be precise, yu p is the displacement of the top vertex point of the beam at the 
free end. 



w 




Vi 




V2 









Hup 




Vpi 




V P 2 









Figure 3.8. Augmented MIMO plant 



Since the system is modeled as a three-input-three-output system, the 
frequency response function (FRF) G(iu) is a 3 x 3 matrix with each element 
Gij(iuj),i,j = 1,2 and 3, corresponding to a particular combination of the 
input and the output; z.e., 

Gij (io;) = G y . Uj = Yi(\u) /U j(\u), 



(3.33) 
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Figure 3.9. Identified model (solid) with measured data (dotted) 



where y 1 = y Up , y 2 = v pl and y 3 = v p2 , and m = w, u 2 = Vi and u 3 = v 2 . 
Yi(iu) and Uj(iu) are the Fourier transforms of yi and Uj respectively. These 
FRFs are determined by applying a sinusoidal chirp of varying frequency 
(from 5 — 250Hz) to the piezoelectric actuators (including the central patch 
corresponding to the disturbance term w) and measuring the corresponding 
output signals yu p , v P i and v p2 . The inputs are generated using a standard 
HP signal generator and the output signals yu p , v P i and v p2 were measured 
using a Polytec laser scanning vibrometer (PSV-300). The vibrometer PSV- 
300 also includes software which computes the ratio of input and output FFTs 
to provide G(iu). 

In Figure 3.9, the FRFs Gij(iu),i,j = 1,2 and 3 are plotted. It is apparent 
from the plots that all the FRFs have three resonance frequencies in the 
plotted frequency region, and the resonance frequencies are more or less the 
same for all the FRFs. 

The goal here is to have good fits of the form (3.22) for the collocated 
FRFs G VplVl and G Vp2V 2 , and good fits of the form (3.2) for all the other 
FRFs. It is known, see [142] and [187], that the beam under consideration has 
a multivariable state space model of the form 
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(3.34) 

(3.35) 

(3.36) 



where 



x(t) = Ax(t) + B w w(t ) + B v V(t ) 

Vtip GyX(t^ 4” Dy W W(t ) + DyyV (t) 

V p (t) = C v x(t) + B vw w(t) + D vv V (£), 
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V(t) = [vi '^ 2 ] / 
Vp{t) = [^pl %>2 ] • 



(3.37) 



(3.38) 



(3.39) 

(3.40) 

(3.41) 



The terms D yw , D vw , D yv and D vv are matrices of appropriate dimensions. 

Using the Laplace transform it can be verified that (3.34)- (3.36) enforces 
the structure (3.22) on the collocated FRFs G VplVl and G Vp2V2 , and the struc- 
ture (3.2) on the other non-collocated FRFs. The model structure (3.34)-(3.36) 
also forces all the FRFs to have the same set of poles. 

Note that the given FRF data has only three resonance frequencies, see 
Figure 3.9. In other words the given FRF data takes into account only the first 
three modes of vibrations. The other higher order modes have been discarded 
(or not taken into account) as they are beyond the frequency region of interest. 
Since the given data includes only three modes, it suffices to set TV = 3, or 
consider a third order model of the form (3.34)-(3.36). A standard method 
to determine the model parameters , /3 k , , u>k , (k} k =i and the feed- 

through terms D yw , D yv , D vw and D vv is to choose them as the minimizers 
of the cost function, 



3 M 



• M = EE 

i,j=l k=l 



Gfj(ivk) ~ Gij(iuj k ) 
Gij (icjp) 



(3.42) 



where G^(icj) are the FRFs corresponding to the multivariable state space 
model (3.34)- (3.36) and denotes the frequency points where the non- 
parametric FRFs Gij( iu) are measured. 
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Minimizing the cost function A4, (3.42), involves using a computationally 
complex non-linear search. Due to the enormity in the number of parame- 
ters to be estimated, the cost function M may have numerous local minimas. 
Therefore a poor initial guess for the parameters , (3k , 7/c , , Ck}^= i, 

D yw , D yv , D vw and D vv (z.e., parameter values which give poor fits for the 
FRF data), the non-linear search may land in a local minima which also gives 
a poor fit for the FRF data. An initialization which in itself gives a decent fit 
to the FRF data would lead the non-linear search to the parameter values that 
give a good fit for the FRF data. One approach to get such initial values is to fit 
the non-parametric data first using subspace methods, [128]. Using subspace 
methods one can obtain a multivariable state-space model, of the same dimen- 
sions as (3.34)- (3.36), for the FRFs Gijfioj). Even though subspace methods 
are black box methods, not conforming to any particular model structure, as 
the cantilever beam is known to possess a model of the form (3.34)-(3.36), 
the subspace fit would also resemble the models (3.34)-(3.36) to a large ex- 
tent. From the subspace fit one can then extract good initial estimates for 
parameters ^ 2 , 7/c, CfcjfcLu D yw, D yv , D vw and D vv . For details 

of subspace methods the reader is referred to [187] or [128]. 

In Figure 3.9 the parametric model of the form (3.34)- (3.36) estimated 
from the FRF data is plotted along with the measured data. It is apparent 
from the plots that the estimated parametric model gives a good fit for the 
non-parametric data. 

3.7.3 PPF Controller Design 

The tip displacement y t i p of the beam gives a measure of the amplitude of 
vibrations in the beam. As w represents disturbance, the FRF G ytipW ( io;) = 
Yap ( iuj) /W ( iu) is a good indicator of the effect of noise on the beam. A well 
damped FRF G ytipW ( io;) would imply a well damped system. Hence, a con- 
troller is designed such that the closed- loop FRF Gci, yw ( i^) corresponding to 
the input w and output y t i p is well damped. Alternatively stated, a controller 
is designed such that the poles of the closed- loop FRF Gci, yw ( i^) are well 
inside the left half plane. 

The PPF controller to be designed is a two-input two-output multi- 
variable controller, which forms a feedback loop connecting the outputs 
V p (t) = [v p i v P 2 \ to the inputs V(t) = [v\ V 2 ]' of the system. Adhering 
to the notations in Section 3.6 the PPF controller to be designed is of the 
form 



x(t) = Ax(t) + rv p (t ) 
V(t) = f'x(t), 



where 



(3.43) 

(3.44) 
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(3.45) 



r 



0 0 

T~t V pl T~\ V p2 

1 1 1 1 



0 0 

T~i v pl t~iV p 2 

1 3 1 3 



(3.46) 



Note that (3.43)-(3.44) is (3.27)-(3.28) written in first order form with N = 3. 

Using (3.43)-(3.44) in (3.34)-(3.36), the closed-loop system can be written 
in the form 



X(t) = A(fi, A, r)X(t) + B w w{t) (3.47) 

ytip(t) = C(r)X(t) + Dy W w(t ), (3.48) 

where X = [x f x']' and 



r = 



~ j^ v p 1 j~i v pi j~\ v pi 

J~\ v p2 j~\Vp2 j~\Vp2 



(3.49) 



Derivations of the expressions for A(f2, A, r),B,C and D yw from (3.43)-(3.44) 
and (3.34)- (3.36) involve standard (and straightforward) calculations and are 
hence omitted. 

The FRF Gci, yw ( i^) corresponding to the closed-loop system (3.47)-(3.48) 
is the one that is to be damped by appropriately choosing 17, A and r. A suf- 
ficiency condition for the stability of the closed-loop (3.47)-(3.48) is provided 
by Lemma 3.6 in Section 3.6. Lemma 3.6 states that for the closed-loop system 
Gci,yw{ h<^) to be stable the LMI condition (3.32) has to be satisfied. Hence 
the PPF controller design problem can be posed as a constrained optimization 
problem: 

min || Gci, yw (uj) || (3.50) 

n,A,r 

subject to the constraint 



!? 2 -VT 0 

-r f & & r 

0 r d - 1 



> 0 , 



(3.51) 



where 




3.7 Experimental Implementation of PPF Control on an Active Structure 



59 



1 ^2 1 ^3 1 " 

2 & 2 2 \ 1 



(3.52) 



Q is the 3x3 diagonal matrix with resonance frequencies of the beam as its 
diagonal elements, and D denotes the feed-through term D vv of the model, 
see (3.36). The design variables T, A and Q are 2 x 3, 3 x 3 and 3x3 matrices 
respectively, with P as defined in (3.49) and the matrices A and Q being 
diagonal. In (3.50) || • || denotes a norm. Here the Hoo-norm is considered as a 
performance measure for the controller design. An H 2 optimal PPF controller 
has also been designed and tested. 

Another approach to design the PPF controller would be to place the 
poles {P£} 4N of Gci, yw (s ) in a desired region of the left half plane. One way 
to design such a controller is to choose (!2, A, T) such that 



4N 4N 

V = aY / \Pk-Re(Pk)\ 2 +P E \Pk-Pk-i \ 2 (3.53) 

k= 1 i= 2, i even 



is minimal under the constraint (3.51). In (3.53), } 47V denotes a set of 

prespecified negative real numbers, Re(-) denotes the real part of a complex 
number, and a and [3 are positive scalars. The first term in the right-hand side 
of the cost function V, (3.53), tries to bring the real part of the closed- loop 
poles {P£} 4N close to the set of prespecfied real values {Pjt} in the left half 
plane. The second term keeps the pair of closed-loop poles that correspond to 
the same open- loop pole close to each other (see Figure 3.10). 
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Figure 3.10. Pole optimisation procedure 



As in the case of identification, both optimization problems (3.50) and 
(3.53) involve using a non-linear search, which needs to be initialized amicably. 
One method to initialize the non-linear search would be to solve the LMI (3.51) 
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using well known LMI solvers like SeDuMi, or the MATLAB® LMI tool box, 
to obtain a feasible solution. Although such an initial condition would be 
valid, its use may result in a controller gain || FT' || which is too large. A 
prohibitively large controller gain may result in implementation difficulties. A 
simple way to obtain initial values is to solve for r by solving the LMI (3.51) 
for certain fixed values of Q and A. Past research (see [177] and [58]) and 
experience with PPF controllers suggest that good PPF designs have their 
resonance frequencies Q close to the system resonance frequencies Q. Taking 
this cue, we let 

Q = kQ (3.54) 

A = —, (3.55) 

K 

where k is a scalar reasonably close to one (for example k = a/ 2), as our initial 
guesses for Q and A. Using the initial values (3.54) and (3.55) a meaningful 
initial guess for F can be obtained in the following fashion: Let and 
denote vectors orthonormal to IF 1 = [^ 1 W r ^ 2 ] f and IF 2 = [iF^ 1 see 

(3.52). Note that due the linear independence of IF 1 and IF 2 (which is assumed) 
there exists C\ and C 2 such that IF 3 = [1F3 1 #3 2 }' is equal to CAF 1 + C 2 lF 2 • 
Let Fq = [aq^/ ol 2 &\ 0] denote a candidate F for the initial value, where 
aq and 0^2 are real valued parameters. It is easy to note that 

aq 0 0 

&'r c = 0 a 2 0 . (3.56) 

C\OL\ C 2^2 0 

Since the feed-thorough term D vv is in general very small, we drop the 
quadratic term r'DT in (3.24) and approximate it by 

Q 2 - &'r c fi- 2 F!j& > 0. (3.57) 

Due to the chosen structure of Fc, (3.57) can be further simplified to 

Q 2 Q 2 > M'M, (3.58) 

where 

OL\ 0 0 

M = 0 ^2 0 . (3.59) 

_CiaifC 2 a 2 f 0_ 

Note that the right-hand side of (3.58) is a constant diagonal matrix, therefore 
choosing aq and a 2 such that the bound (3.58) is satisfied would give an initial 
guess for Fc- However, it must be stressed that the approximation introduced 
in (3.57) could fail for large values of aq and a 2 . Hence, it is worth checking 
if the LMI (3.24) is satisfied for the chosen initial guesses. 
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Figure 3.11. Schematics of the experimental setup 



Results of Pole Placement Controller Design 

In this subsection, the results pertaining to the effectiveness of the pole opti- 
mized PPF controller are presented. In Figure 3.12, the two-input two-output 
multivariable pole optimized PPF controller is plotted. Note from the plots 
that the FRFs of the PPF controller have a good roll-off in the high frequency 
regions, and they also do not have large gains. These twin properties of the 
PPF make them easily implement able. By construction, see (3.21), the di- 
agonal entries of the PPF controller correspond to collocated FRFs and the 
cross diagonal terms refer to the non-collocated FRFs. It can be noted from 
Figure 3.12 the non-collocated FRFs are identical, which is expected due to 
(3.21). 

Schematics of the experimental setup used to implement the pole opti- 
mized PPF controller, and all subsequent controllers, is illustrated in Fig- 
ure 3.11. The controller was designed, and simulations were performed, in 
MATLAB® Simulink®. The controller was then downloaded onto a dSPACE 
DS-1103 rapid prototyping system. Sampling frequency of the DSP system 
was set to 20 KHz, and low-pass antialiasing and reconstruction filters with 
cut-off frequencies of 10 KHz were added to the system as shown in Fig- 
ure 3.11. Voltages induced in the two piezoelectric sensors were measured 
through high-input-impedance buffer circuits, to reduce the low-frequency dis- 
tortions caused by the finite impedance of the measurement device. 

To determine a controller, the desired real parts of closed-loop poles in 
(3.53) were placed at —8, —66 and —120 respectively. The real parts of the 
actual closed-loop poles were found to be —9.7494, —60.8264 and —105.9135 
respectively. 

To evaluate the damping introduced in the system due to the PPF con- 
troller, both the open- loop and closed- loop FRFs G ytipW ( io;) and Gci,y tipW {^) 
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Frequency (Hz) 



Figure 3.12. PPF Controller using the pole optimization 



are plotted in Figure 3.13 and 3.14. In Figure 3.13 simulated values of 
Gy Up w(iu) and Gci, ytipW {^) are plotted, i.e., G ytipW (iuj ) obtained from the 
model (3.34)-(3.36) and Gci, ytipW ( i^) obtained from (3.47) and (3.48) with 
the parameters 12, A and r set to their optimal values. In Figure 3.14 experi- 
mentally determined G ytipW (iw) and Gci, ytipW {^) are plotted. Experimentally 
determined FRFs refer to the case where the piezoelectric patch correspond- 
ing disturbance input w in the cantilever beam setup is driven by a chirp 
signal and the corresponding output yu p is recorded for both the closed- and 
the open-loop systems. The FRFs G ytipW ( icj) and Gci, ytip w( i^) are then cal- 
culated from the recorded input-output data, as done in the non-parametric 
identification. Simulations and experimental results are fairly similar, thus val- 
idating the use of the model (3.34)- (3.36) for the beam. The plots also suggest 
a satisfactory damping in the magnitude of the closed- loop FRF Gci, ytipW ( i^) 
at the resonance frequencies. 

In Figure 3.15 the open- loop poles and the closed-loop poles are plotted. It 
is evident that Gci, ytipW { i^) is much better damped compared to G ytipW (iu). 

In practice, due to wear and tear and also due to other factors such as 
surrounding temperature, etc. the material properties of the structure and 
the piezoelectric patches tend to change. This leads to shifts in the resonance 
frequencies of the beam. A good controller design must be robust enough to 
provide good damping even under changed circumstances. In order to check 
the robustness of the PPF design artificial shifts in the resonance frequencies 
of the beam were brought about by adding an extra mass at the free end of 
the beam. In Figure 3.16 the measured open- loop and the closed- loop FRFs 
Gyupwfa) and Gci, ytip w( i^) for the loaded system are plotted. It must be 
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13. Simulated open- and closed- loop frequency response using the pole 




Figure 3.14. Measured open- and closed- loop frequency response using the pole 
optimization 
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Figure 3.15. Open- and closed- loop pole map for the pole placement optimization 
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stressed that, the PPF controller used in the closed- loop system is the same 
PPF controller that was used for the unloaded Cantilever beam. The plots 
suggest that the PPF controller designed using pole optimization is robust 
with respect to shifts in the resonance frequencies. 

Finally to illustrate controller performance in the time-domain, a pulse- 
shaped disturbance voltage was applied to the disturbance piezoelectric patch. 
The resulting tip displacements in open- and closed-loop were recorded for the 
beam with and without the tip mass. These time-domain results are plotted 
in Figure 3.17 and clearly demonstrate the effectiveness of the controller in 
terms of disturbance rejection and its relative immunity to uncertainty in the 
structural dynamics of the system. 

iToo PPF Controller 

Finally, this section is concluded with the results of using an optimized 
PPF controller. FRF of the optimized PPF controller is plotted in Fig- 
ure 3.18. As in the previous case the PPF controller designed by minimizing 
the Hqq norm also has a good roll-off in the high frequency regions. In order 
to limit the level of control signals applied to the two piezoelectric actuators, 
the exogenous output yu p is augmented with weighted versions of v\ and V 2 . 
The weights are chosen carefully to trade off between the achievable damping 
and the necessary control signals. This is a standard approach [206], hence no 
more details are included here. 
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Figure 3.16. Measured open- and closed- loop frequency responses using the pole 
optimization, with and without the mass 






Figure 3.17. Measured disturbance response using pole optimization, with and 
without the mass 
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Figure 3.18. PPF Controller using the Hoo optimization 



In Figures 3.19 and 3.20 the closed-loop and the open-loop FRFs are plot- 
ted. The plots suggest a good damping of the closed-loop due to the PPF 
controller. The plots also suggest good agreement between the simulated and 
experimental FRFs. 

In Figure 3.21 the poles of the open- and the closed- loop systems are 
plotted. As in the earlier case the plot suggests a good deal of damping in the 
closed-loop system due to the PPF controller. 

Robustness of the optimized PPF controller is tested by using it on 
the loaded cantilever beam. The plots 3.22 suggest that the minimized 
PPF controller is also robust to changes in the resonance frequencies of the 
beam. 

Finally, to observe the performance of the optimized PPF in the time- 
domain, similar time-domain tests are performed on the system. The results 
are illustrated in Figure 3.23. 



3.8 Self-sensing Techniques 

In a typical active vibration control application, piezoelectric elements are of- 
ten used as actuators, or sensors. In this case, the piezoelectric device performs 
a single function; either sensing, or actuation. The piezoelectric self-sensing 
actuator, or sensori- actuator, on the other hand, is a piezoelectric transducer 
used simultaneously as a sensor and an actuator. This technique was devel- 
oped concurrently by Dosch, Inman, and Garcia [54] ; and Anderson, Hagood 
and Goodlife [9] , who made the observation that with the capacitance of the 
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Figure 3.23. Measured disturbance responses using the optimization, with and 
without the mass 



piezoelectric device known, one can simply apply the same voltage across an 
“identical” capacitor and subtract the electrical response from that of the 
sensori-actuator to resolve the mechanical response of the structure. 

The key idea, here, is to replace the function of a sensor in the feedback 
loop by estimating the voltage induced inside the piezoelectric transducer, 
v p . Since this voltage is proportional to the mechanical strain in the base 
structure, the estimated signal would provide a meaningful measurement for 
a feedback compensator. Furthermore, by estimating v p , one would effectively 
replace the role of the collocated piezoelectric transducer in Figure 3.1 by 
the additional electronic circuitry. In this way one would expect to design 
feedback controllers that possess appealing properties associated with com- 
pensated collocated systems, such as those discussed in previous sections. 

Two realizations for the piezoelectric sensori-actuator, as proposed in [9], 
are sketched in Figures 3.24 and 3.25. The two circuits have rather similar 
functions; they use a signal proportional to the electrical charge or current 
and subtract that from the signal proportional to the total charge or current 
to produce a signal proportional to the mechanical strain, or its derivative. 
This signal is then used for feedback. 

In the strain measurement circuit of Figure 3.24, assuming that the leakage 
resistors R\ and R 2 are very large and that the gain of each op-amp voltage 
follower is unity, we may write: 
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Figure 3.24. The piezoelectric-based sensori-actuator, generating an estimate of 
the mechanical strain 
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where v p and C p are respectively, the voltage induced in and the capacitance 
of the piezoelectric transducer (refer to Figure 2.16 (b)). The voltage v p is 
proportional to the mechanical strain. Subtracting v 2 from v\, we obtain: 



V s = Vi - v 2 

( C r 

VCi + a 



Cp A 

c 2 + c p ) 



V + 



Cp 

c 2 + c p 



v p . 



Hence, if Ci = C 2 and C r = C p , equation (3.60) reduces to 



s 



Cp 

C‘2 + Cp 1 



(3.60) 



v 



(3.61) 
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Figure 3.25. The piezoelectric-based sensori- actuator, generating an estimate of 
the strain rate 



Therefore, under the above “ideal assumptions” , the estimated voltage is 
proportional to v p . Now, consider the sensori- actuator in Figure 3.25. The 
voltage v is applied to both the piezoelectric transducer and the reference 
capacitor, C r . A current, i e , flows through the upper path, while i p flows 
through the lower path. Each signal is converted to a voltage using an op- 
amp. The two signals are differenced, resulting in a voltage proportional to 
the derivative of v p ; z.e., the strain rate. To be more precise, if the two resistors 
Ri and R 2 are both equivalent to i7, v s is found to be: 

v s = R(C r — C p )sv + RC p sv p . 

Again, it can be observed that if C r = Cp, the first term will disappear, 
and v s will be proportional to the strain rate. For practical reasons, however, 
very often the capacitive and resistive elements are chosen differently; see [9] 
and [54] for more details. 

The two sensori- actuator schemes in Figures 3.24 and 3.25 should perform 
well under ideal assumptions. Having estimated or perhaps v p , the signal 
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produced by the sensori-actuator can now be used for feedback. Several appli- 
cations for this method have appeared throughout literature (see for example 
[81, 200, 106, 21, 185, 4]. In practice, however, there are a number factors that 
limit the performance of the sensori-actuator, the foremost being the choice of 
the reference capacitor C r that is directly related to the size of piezoelectric 
capacitance C p . The piezoelectric properties are influenced by variations in 
environmental conditions and operation. This requires a continual effort to 
tune the circuits in Figures 3.24 and 3.25. 

A primary obstacle for implementation of the piezoelectric sensori-actuator 
is the difficulty in obtaining an accurate estimation of the capacitance of the 
piezoelectric device, C p . This may not severely affect the open- loop perfor- 
mance of the sensori-actuator, however, if v s is used as measurement for feed- 
back, such variations may destabilize the closed-loop system 2 . An attempt to 
address this problem was made in [38, 191, 36, 1], where the authors suggest 
an adaptive sensori-actuator implementation based on the LMS algorithm 
[195, 55]. 

The sensori-actuator is a linear estimator that generates an estimate of 
the strain signal, or its derivative. The structure of the estimator is rather 
crude and largely dependent on the added electronic circuitry. Often a nominal 
model for the underlying system is at hand. Therefore, it should be possible to 
construct better estimates of the required signals using an optimal estimation 
method such as a Kalman Filter [114, 7]. The issue of uncertainty associated 
with the varying piezoelectric capacitance can then be addressed using the 
recent advances in robust state estimation and Kalman Filtering (see [156] 
and references therein). 



2 In fact, it can be shown that the transfer function estimated by the self-sensing 
circuit is G vv (s) + <5, where S is proportional to C p — C r . This additional feed- 
through term does not alter poles of the open-loop system. However, it does 
perturb the open-loop zeros, and this could be detrimental to the closed-loop 
performance and stability of the system. 
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Piezoelectric Shunt Damping 



4.1 Introduction 

Piezoelectric shunt damping is a popular technique for vibration suppression 
in smart structures. As illustrated in Figure 4.1, techniques encompassed in 
this broad description are characterized by the connection of an electrical 
impedance to a structurally bonded piezoelectric transducer. Such methods 
do not require an external sensor, may guarantee stability of the shunted 
system and do not require parametric models for design purposes. 

The first goal of this chapter is to provide a brief review of the various 
techniques categorized as piezoelectric shunt damping. One particular sub- 





Figure 4.1. A piezoelectric laminate structure disturbed by an external force w. 
The resulting vibration d is suppressed by the presence of a shunt impedance Z(s). 
The poling direction of the transducer is indicated by the shaded arrows. 
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category, the so-called resonant shunts, has been the focus of significant de- 
velopment over the past decade. Resonant shunt damping circuits, comprised 
of inductors, capacitors, and resistors, are simple to design and can signifi- 
cantly augment the damping of lightly-damped flexible structures. The fore- 
most difficulty associated with resonant shunt circuits is the requirement for 
impractically large inductance values. In Section 4.4 the synthetic impedance 
is introduced as a simplified implementation technique. 

In the next chapter, the structural influence of a shunted piezoelectric 
transducer is studied and expressions revealing the dynamics of a piezoelectric 
laminate structure are derived. The technique of piezoelectric shunt damping 
is shown to be equivalent to collocated strain feedback control. A simple re- 
lationship is shown to exist between a connected shunt impedance and an 
equivalent strain feedback controller, and vice versa. 



4.2 Passive Shunt Damping 

In the following subsections, the various passive shunt damping techniques 
depicted in Figure 4.2 are discussed in detail. 

4.2.1 Passivity Defined 

An electrical shunt impedance is said to be passive if and only if it does not 
supply power to the system. In mathematical terms, this means that [6] 

oo 

v(t) • i(t) > 0, (4.1) 

where v(t) and i(t) are the voltage and current as illustrated in Figure 4.1. 
For linear systems, passivity of the impedance Z(ju) is guaranteed if [6] 

$t(y(juj) • i*(ju ) ) > 0 or $t(Z(jw)) > 0 Vo;, (4.2) 

where i*(jw) is the complex conjugate of 

The foremost benefit of passivity is that stability of the shunted system 
is guaranteed. A stability analysis of passively shunted piezoelectric laminate 
structures is presented in Section 5.5. 

In the context of shunt damping and vibration control, the reader should 
be aware that the definition of electrical passivity can vary between authors. 
In some discussions, passivity not only implies that the above conditions hold, 
but also that the impedance is constructed solely from physical components 
with no external power supply. This meaning is occasionally used by some au- 
thors while discussing the characteristics of particular shunt damping circuits; 
for example, a shunt circuit that operates autonomously without a power sup- 
ply is commonly referred to as “passive” . In this book we adhere to the stan- 
dard definition as discussed in [6]. For example, in Section 4.4 op-amp circuits 
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Figure 4.2. Passive piezoelectric shunt damping techniques discussed in Section 
4.2 



are introduced as a means for implementing impractically large inductance 
values. Although the impedance is synthesized using active components, the 
terminal voltage and current satisfy the conditions in Equations (4.1) and 
(4.2), thus it is referred to as “passive”. 
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Figure 4.3. A piezoelectric transducer with an RL shunt and the equivalent elec- 
trical circuit 



4.2.2 Linear Techniques 

For the purposes of this discussion, a linear shunt circuit is defined as any 
impedance with a linear current-to- voltage relationship over the bandwidth of 
interest. By this definition, slowly varying impedances such as adaptive shunts 
(Chapter 8) and shunts implemented by switched mode amplifier (Section 6.7) 
are classed as linear. 

In reference [82] a simple resistance was applied to a structurally bonded 
piezoelectric transducer in order to damp vibration. Its effect on structural 
dynamics was shown to be equivalent to that of a light visco-elastic damping 
treatment [82]. Although simple, resistors are rarely used as they offer only a 
small amount of damping; typically only a few dB when applied to a lightly 
damped Aluminum structure. Future transducers utilizing high dss electro- 
mechanical coupling factors may provide better performance. 

Capacitive shunting [49] is another simple but poorly performing tech- 
nique. Adding capacitance to the terminals of a piezoelectric transducer varies 
its effective stiffness. Although capacitive shunts have not proven useful in 
vibration control, they may find application where a slight adjustment in res- 
onance frequency is required. Capacitive shunts and their equivalent strain- 
feedback interpretation are discussed further in Section 5.3. 



4.2.3 Resonant Single-mode Shunt Circuits 

The first shunt circuit of any kind should be credited to Forward [68], who 
in 1979 proposed the idea of inductive (LC) shunting for narrow band reduc- 
tion of resonant mechanical response. In particular, he demonstrated that the 
effect of inductive shunting was to cancel the inherent capacitive reactance 
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Figure 4.4. Circuit implementation of (a) a grounded and (b) a floating inductor 



of a piezoelectric transducer. Later, Hagood and von Flotow [82] interpreted 
the operation of a resonant shunted piezoelectric transducer as analogous to 
a tuned mass damper, in which a relatively small second order system is ap- 
pended to the dynamics of a larger system. Moreover, they addressed the 
situation in which a resistive element is added to the shunt network, result- 
ing in an RLC tuned circuit. This system, and its electrical equivalent are 
depicted in Figure 4.3. The resulting RLC circuit is tuned to a specific reso- 
nance frequency of the composite system. That is, if the vibration associated 
with the £ th mode is to be reduced, then L is chosen as: 



C p m 

By adopting a proper value for R, the resonant response at and in the 
vicinity of uj£ can be reduced. However, we should keep in mind that due to 
the passive nature of the shunt there will be hard constraints on the level of 
performance that this approach can achieve. 

A number of techniques can be used to determine an optimal value for the 
resistive component of the shunt impedance. One possibility is to set up an 
optimization problem to minimize the H 2 norm of a transfer function associ- 
ated with the structure as a function of the parameter R. Such a procedure is 
employed in reference [16] and if performed properly, will result in the place- 
ment of poles so that structural damping is augmented. Alternatively, one 
could employ a pole placement method, such as explained in reference [82]. 

A parallel circuit variation where the resistance and inductance are placed 
in parallel rather than series was proposed by Wu [197]. Although the two 
circuits, series and parallel, achieve similar performance, it is claimed that 
the parallel structure is less sensitive to sub-optimal resistance values. 
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The foremost difficulty in implementing resonant shunt circuits is the need 
for very large inductive elements. As an example, to implement a single- mode 
RL shunt for control of a vibratory mode of 11 Hz using a piezoelectric trans- 
ducer with capacitance of C p = 101 nF requires an inductor of 2072 H. Such 
an inductor can only be synthesized using active electronic circuitry. Two 
possible implementations are demonstrated in Figure 4.4 1 . 

As illustrated in Figure 4.4 (a), a grounded inductor can be synthesized 
by two op- amps. The equivalent impedance observed at the terminals of this 
Gyrator is 

v _ Z1Z3Z5 

^ in — r? r? ' 

■^2^4 

Allowing Z 4 to be a capacitor, C, and replacing other impedances with resis- 
tors, the resulting circuit will act as an inductor L = kC , where 

_ R1R3R5 

~ R‘2 ' 

This circuit is all that is needed for implementation of single-mode RL shunts. 
To implement certain classes of multi-mode shunts, however, one would need 
to implement floating inductors. This can be done using four op-amps, shown 
in Figure 4.4 (b). 

Replacing the inductor needed in a single-mode RL shunt by a Gyrator is 
practical. However, particular care must be taken in implementing such cir- 
cuits. Voltages generated across piezoelectric transducers could be very large, 
particularly when disturbances acting on the structure have frequency content 
that is very close to resonance frequencies of the base structure. Subsequently, 
high-voltage op-amps are to be used for this purpose. The designer should also 
be careful to avoid saturation of internal nodes of the system. For a single- 
mode RL shunt this is straightforward, due to the simplicity of the required 
circuit. For multi- mode circuits such implementation becomes overly complex 
and impractical. This issue will be further discussed in Section 4.4. 

4.2.4 Resonant Multi-mode Shunt Circuits 

Extension of the single-mode shunt damping technique to allow for multiple- 
mode vibration suppression has been the subject of intense research in recent 
years. A trivial solution is to attach a number of piezoelectric transducers 
to a structure, each one shunted by an RL circuit tuned to a specific mode. 
This is clearly not a viable option as there is only a limited amount of surface 
area over which the transducers can be mounted. The main focus in this area, 
therefore, has been on devising multiple-mode vibration damping methods 
which use a single piezoelectric transducer. The remainder of this chapter is 
devoted to the study of this issue. 

1 It is possible to reduce the size of necessary inductors by adding an extra ca- 
pacitance in parallel to the piezoelectric transducer. This issue will be further 
elaborated on in the next chapter. 
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Figure 4.5. Two- mode shunt damping circuit [198] 



Current-blocking 

The first multi-mode techniques were proposed by Wu and his co-authors 
[199, 198, 203]. Their idea is centered around using an RL (either parallel, or 
series) shunt for each individual mode, and then inserting “current blocking” 
LC circuits into each branch. The electric shunt circuit for a two-mode system 
is depicted in Figure 4.5. 

If vibration of the first two modes of the base structure are to be reduced, 
then L\C V is tuned to uj% while L 2 C P is tuned to CJ 2 - Furthermore, L\C\ 
is tuned to CJ 2 while L 2 C 2 is tuned to u)\. Therefore, L\Ri and L 2 R 2 are 
effectively separated at and CJ 2 - For three modes, two current blocking 
circuits are inserted inside each branch, and so on. 

The main difficulty with this method is that the circuit order increases 
very rapidly with the number of modes to be shunt damped simultaneously. 
This seriously complicates the task of implementing the required circuits. 




Figure 4.6. Simplified two- mode shunt 
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Figure 4.7. Hollkamp circuit 



The number of components needed in this multi-mode shunt circuit can 
be almost halved by using straightforward circuit theory. For example, the 
two- mode shunt in Figure 4.5 can be reduced to that in Figure 4.6. For this 
circuit to work properly, its components must be chosen according to 



U 

L 2 



1 

(LiL 2 + L 2 L 2 ~ L\L 2 — {jj\L\L 2 Ij 2 C 2 
(Li ~ L 2 ) (l - cc ; 2 2 L 2 C 2 ) 



(4.3) 

(4.4) 



and 



> bll 

to 

II 

_ & ^ 


(4.5) 


l 2 = 

ulc 2 


(4.6) 



Here, C 2 is an arbitrary capacitor used in the current blocking network. 

This class of multi-mode impedances can be tuned to damp a small number 
of modes in an effective way. An example of this is contained in [16] where the 
second and third modes of a simply supported beam are reduced in magnitude 
by 18.9 dB and 19.1 dB. 

Hollkamp 

Hollkamp [89] suggested a specific resonant shunt structure, depicted in Fig- 
ure 4.7. The shunt circuit consists of a number of parallel RLC shunts, with 
the very first branch being an RL circuit. For one mode, Hollkamp’s circuit 
reduces to the one proposed by Hagood and von Flotow. However, for each 
additional mode, an RLC branch has to be added. When an extra branch 
is added, the previous resistive and inductive elements must be re-tuned to 
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Figure 4.8. Current- flowing multi-mode shunt circuit [17] 



ensure satisfactory performance. No closed form tuning solution has been pro- 
posed for this technique. However, in [89] values for the shunt circuit electrical 
elements are determined by numerical optimization of a non-linear objective 
function. Given that all circuit elements are to be determined numerically, for 
a large number of modes this procedure may result in a complicated optimiza- 
tion problem. Nevertheless, this method has been applied to a cantilevered 
beam in [89] , in which vibration of the second and third modes were reduced 
by 19 and 12 dB respectively. 



Current-flowing 



The current- flowing shunt circuit [17] was introduced as a means for sim- 
plifying the implementation of high-order multi-mode shunt circuits. Shown 
in Figure 4.8, a current-flowing shunt circuit requires one parallel branch for 
each structural mode to be controlled. The current-flowing network LiCi in 
each branch is tuned to approximate a short circuit at the target resonance 
frequency whilst approximating an open circuit at the frequencies of adjacent 
branches i.e. 



Li 



1 

u&i 



• • • , l n 




(4.7) 



The remaining inductor and resistor in each branch LiRi, are tuned to damp 
the i th target structural mode in a manner analogous to single-mode shunt 
damping i.e., 



Lt 



1 



Ln = 



00 



N^P 



(4.8) 



The two inductors in each branch can be combined, resulting in a series 
RLC circuit in each parallel arm of the shunt circuit. This further simplifies 
the task of implementing the shunt circuit. The total shunt admittance is then 
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Figure 4.9. The series-parallel shunt damping circuit of [17]. 

N 

Y(s) = J2Yi(s) 

i = 1 

= (V Lj)s 

S s* + (Ri/LJs + (l/LiCiY 

where 

Li ~ Li T Li . 

Compared to the circuit proposed by Wu and his co-authors [199, 198, 203] , 
the resulting shunt circuit is of a considerably lower order. Furthermore, in 
comparison with the technique proposed in [89], this is a more systematic 
way of designing a shunt impedance circuit. The only parameters left to be 
chosen are the resistors i?i, . . . , Rn- These parameters can be determined via 
optimization or tuned experimentally. 

Series-parallel 

The electrical dual of a current-flowing circuit, the so-called series — parallel 
circuit shown in Figure 4.9 was proposed as a method for reducing inductive 
component values in [67]. 

Each network in series CiLiLiRi contains two sub- networks: a current- 
blocking network CiL^ and a parallel single-mode damping network LiRi. 
Both the current-blocking and damping networks, CiLi and LiRi, are tuned 
to the same target resonance frequency u At this frequency, the current- 
blocking network has an extremely large impedance. All of the remaining 
current-blocking networks, tuned to other structural resonance frequencies, 
have a low impedance at Ui. A voltage applied at the terminals results in a 
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current that flows freely through the detuned current-blocking networks but 
is forced to flow through the active damping network. In this way, the circuit 
is decoupled so that each damping network LiRi can be tuned individually to 
a target resonance frequency. As described above, the series-parallel structure 
offers no great advantage over comparable techniques. Benefits arise from a 
suitable choice in the arbitrary capacitances C{. The recommended capaci- 
tance value is 10 — 20 times that of the piezoelectric capacitance. In this case, 
the current blocking inductors become significantly smaller than the damp- 
ing inductors. When the circuit is simplified by combining the damping and 
current-blocking inductors in each series network, the resulting single inductor 
is a fraction of that required in other single- or multi- mode circuits. 

Control Orientated 

All of the multi-mode techniques discussed thus far are more or less direct 
extensions of the original single- mode circuits. A new approach to the design of 
passive piezoelectric shunt damping circuits was presented in [140]. By viewing 
the electrical impedance as parameterizing an equivalent collocated strain 
feedback controller, a shunt impedance can be found by working backwards 
from an effective feedback controller. Under certain conditions, the passivity, 
and hence stability of the shunted system can be guaranteed [140]. Present 
controller designs have benefits similar to that of current-flowing circuits, they 
are low in order, easy to tune, and suitable for modally dense systems. The 
design and analysis of these shunt circuits is discussed in Section 5.5. 



4.2.5 Non-linear Techniques 

In an attempt to eliminate the need for large inductors, literature has also 
developed on the so-called switched shunt or switched stiffness techniques 
[41]. Three major subclasses exist where the piezoelectric element is switched 
into a shunt circuit comprising either of: another capacitor [49], a resistor [37], 
or an inductor [163]. The required inductance is typically one-tenth of that 
required to implement a simple LR resonant shunt circuit designed to damp 
the same mode. Similar to virtual circuit implementation, an external power 
source is required for the gate drive and timing electronics. 

In reference [148] a hybrid systems approach is taken to derive optimal 
switching rules for the aforementioned switching shunts. It is shown theoreti- 
cally and experimentally that the empirical switching laws can be improved. 
Although the optimal solution is numerically intensive to compute in real- 
time, a simplified quasi-optimal solution is easily implemented. 

The Strain Amplitude Minimization Patch (STAMP) damper was pre- 
sented as a self-powered shunt vibration controller in references [107] and [108]. 
Achievable performance is less than switched inductor designs but without the 
requirement for a power supply. 
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Figure 4.10. Active piezoelectric shunt damping techniques 



A better performing autonomous shunt using two transducers was pre- 
sented in [148]. While one transducer is connected to a switching circuit closely 
resembling a switched inductor network, the second transducer is utilized to 
provide energy for MOSFET gate drives and to correctly time the switching 
events. 

Another non-linear technique, the time varying resistance, was suggested 
for maximizing the loss factor of a shunted system [193]. The foremost draw- 
backs are the complicated design process and complexity of implementation. 



4.3 Active Shunt Damping 

Active shunt impedances cannot be realized using passive physical compo- 
nents. Although passivity, and hence stability is not guaranteed, active shunts 
have been known to provide greater vibration suppression than passive cir- 
cuits. 

The negative capacitor shunt circuit [201, 15] is a simple technique for 
broad-band structural damping. By treating the internal piezoelectric voltage 
source as a supply, and the shunt impedance as a load, the traditional concept 
of maximum power transfer can be applied. The optimal impedance is equal 
in magnitude to the source impedance, but opposite in phase; hence the neg- 
ative capacitor. Detrimentally, the optimal impedance applies large control 
voltages from DC to out-of-bandwidth frequencies. An effective technique for 
stabilizing and toning down the control effort of a negative capacitor is pre- 
sented in [15]. Although negative capacitor shunts are somewhat immune to 
variation in the structural resonance frequencies, any variation in the trans- 
ducer impedance can heavily degrade performance and lead to instability. Due 
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to their simplicity, negative capacitors are discussed further and evaluated ex- 
perimentally in Chapter 9. 

Optimal shunts are discussed at length in Chapter 10. By viewing the 
problem of shunt impedance design as a standard feedback problem. Synthesis 
techniques such as LQG , # 2 , and are readily employed to design opti- 
mal impedances. Optimal shunts can provide excellent vibration suppression 
(greater than 30dB, see Chapter 10) and can be made less sensitive to varia- 
tion in resonance frequencies. The main disadvantage associated with optimal 
shunts is the requirement for a dynamic system model and the dependency 
on transducer capacitance. 

The state-switched active shunt (SSDV) [157] was proposed to increase the 
performance of switched inductor techniques. Rather than switching a series 
inductor-resistor network between open- and closed-circuit states, the network 
is switched between two voltage sources. Damping performance is increased 
but at the expense of additional power requirements. 

The final vibration control technique falling loosely into the category of 
piezoelectric shunt damping is hybrid active-passive control [182]. This tech- 
nique comprises a sensor and piezoelectric actuator driven by a voltage source 
in series with a resonant shunt circuit. There are a number of advantages in 
augmenting an active feedback system with a passive shunt: 1) Damping is 
still provided by the shunt circuit if the sensor fails or if the voltage source 
fails in short-circuit; 2) The system transfer function seen by the controller 
consists of the true structural system with damping augmented by the pas- 
sive shunt. This can simplify control design and provide greater immunity to 
changes in resonance frequency. 



4.4 Implementation of Resonant Shunt Circuits 

For typical values of piezoelectric capacitance and structural resonance fre- 
quency, resonant shunt circuits require inductance values in the order of tens 
or hundreds of Henrys. Inductors of this size must be implemented by virtual 
circuits or otherwise. 



4.4.1 Virtual Circuits 

Virtual inductors and Riordan gyrators [164] have been employed to imple- 
ment large inductance values. They are most useful for the implementation of 
single-mode shunt circuits. In reference [152], a practical single-mode circuit 
is described with the option for online analog tuning. 

The foremost difficulties associated with virtual circuit implementation are 
discussed below. 

• Virtual inductors or Riordan Gyrators are difficult to tune and are sen- 
sitive to component age, temperature, and non-ideal characteristics. The 
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complexity of the required circuitry is illustrated in Figure 4.11 where a 
simplified dual-mode current-blocking network is implemented by virtual 
inductors and Riordan Gyrators. 

• As piezoelectric patches are capable of developing hundreds of volts from 
moderate structural excitations, the entire circuit must be constructed 
from high-voltage components. In the likely case of voltages greater than 
±45 V, the components may be prohibitively expensive. Virtual circuits 
may also contain internal voltage gains. In this case, the range of operation 
is dictated by the greatest circuit node voltage. 

• The minimum number of op-amps required to implement a piezoelectric 
shunt damping circuit increases rapidly with the number of modes to be 
damped. At least 30 op-amps are required to implement a three-mode se- 
ries configuration multi-mode shunt damping circuit with current-blocking 
networks in every branch. In general for this configuration, 2n ± 4n(n — 1) 
op-amps are required to damp n modes. If 5 modes are to be controlled, 
90 op- amps are required. The advent of current-flowing shunt circuits has 
significantly alleviated circuit complexity issues since they require only 2 n 
op- amps to control n modes. 

Despite the associated practical problems, the majority of work within the 
field has adopted virtual circuits for shunt impedance implementation. This 
may be a cause for the slow transfer of piezoelectric shunt damping technology 
into the industrial domain. 

4.4.2 Synthetic Admittance 

From the discussion in the previous section it should be clear that a key 
difficulty associated with passive piezoelectric shunt damping is the need for 
large inductive elements. As demonstrated in Figure 4.11, the necessary circuit 
can be implemented using virtual inductors and gyrators. The complexity of 
such implementation precludes the extension to high-voltage or high-order 
shunt circuits. Here, the synthetic admittance is introduced as an alternative 
that is not limited by shunt circuit complexity or voltage. 

As originally proposed in [59] and [60], the concept is illustrated in Figure 
4.12. The synthetic admittance is a two terminal device consisting of a voltage 
buffer, signal filter V(s), and a voltage controlled current source (VCCS). 
Referring to Figure 4.12 (b), the applied current i z is determined by the 
measured voltage v z . By fixing i z as the output of a linear transfer function 
V(s), the controlled current source can be made to synthesize an arbitrary 
terminal impedance Z(s). In the Laplace domain, 

i z (s) = Y(s)v z (s), (4.9) 

where Z(s) = y— is the desired terminal impedance. 

In some circumstances it is advantageous to utilize a voltage source rather 
than a current source. In this case, referred to as impedance synthesis , Fig- 
ure 4.13 (b), the applied voltage v z is determined by the measured current i z . 
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Figure 4.12. A synthetic admittance circuit (b) employed to implement a shunt 
impedance Z(s). 



By fixing v z as the output of a linear transfer function Z(s), the controlled 
current source can be made to synthesize an arbitrary terminal impedance 
Z(s). In the Laplace domain, 



v z (s) = Z(s)i z (s), (4.10) 

where Z(s) is the desired terminal admittance. 

The choice of configuration, either synthetic impedance or admittance, will 
depend on the relative order of the desired impedance. As implementation of 
improper transfer functions is impractical [101], the choice should be made 
so that the required transfer function Z(s) or Y(s) is at least proper [101]. 





Figure 4.13. A synthetic impedance circuit (b) employed to implement a shunt 
impedance Z(s). 
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Figure 4.14. The experimental setup. An artificial disturbance is introduced by 
applying a voltage v a to a disturbance transducer. The resulting vibration is sup- 
pressed by the presence of a shunt circuit implemented by a synthetic admittance 



The former method, however, is believed to be more advantageous. It is often 
more straightforward to obtain high-precision voltage measurements from a 
piezoelectric transducer, while the current can be supplied with the required 
precision. Another justification for this observation is related to the hysteretic 
behavior of piezoelectric transducers at higher drives, when driven by a volt- 
age source. It is known that a piezoelectric transducer displays negligible hys- 
teretic nonlinearities if it is driven by a current source instead [78, 147, 39]. 

Further discussion on the implementation of synthetic admittance circuits, 
including current and charge amplifiers for piezoelectric loads can be found in 
Chapter 6. 



4.5 Experimental Demonstration 

In this section, a dual-mode current-blocking circuit is implemented by syn- 
thetic admittance to damp the second and third modes of a simply supported 
beam. The resonance frequencies of the experimental apparatus, shown in Fig- 
ure 4.14, are 74.5 Hz and 171.3 Hz respectively. In this experiment two iden- 
tical collocated transducers are utilized, one to introduce an artificial strain 
disturbance, and another to control resulting vibration. Both of the transduc- 
ers are made from PIC 151 piezoelectric ceramic and have a capacitance of 
101 nF. Vibration of the beam at a point on its surface is measured using a 
PSV300 laser vibrometer. 

The inductance values Li, L 2 and L 2 of the circuit shown in Figure 4.6 
were determined from Equations (4.3) to (4.6). The capacitance C 2 , for the 





90 



4 Piezoelectric Shunt Damping 




Figure 4.15. Simulink® diagram of the two- mode shunt 



“blocking circuit” was chosen as lOOnF. The inductor parameters were found 
to be L\ = 43 H, L 2 = 20.9 iJ, and L 2 = 45.2 H. 

To determine optimal values for the resistive elements, a state-space model 
of the composite system was derived and parameterized in terms of R\ and 
R 2 . An optimization problem was then set up to minimize the H 2 norm of the 
transfer function from the applied disturbance voltage to the displacement 
on the surface of the beam. A minimum was found at R\ = 262.75 kfi and 
R* 2 = 550.73/cC. 

The shunt was then implemented using a simple synthetic admittance 
circuit. The synthetic admittance circuit used in this experiment is depicted 
in Figure 6.7 and fully explained in Section 6.6. The admittance transfer 
function Y(s) was implemented digitally on a dSPACE 1103 Digital Signal 
Processor. A block diagram of the required admittance transfer function is 
shown in Figure 4.15. C-code for the DSP was generated automatically using 
the Real Time Workshop for MATLAB®. 

Frequency-domain results are demonstrated in Figure 4.16. The exper- 
imental resonant amplitudes for the 2nd and 3rd modes were successfully 
reduced by 22 and 1 8dB respectively. To demonstrate the effect of additional 
damping on the transient response of the structure, a second set of exper- 
imental data was collected. The experiment consisted of exciting the beam 
structure using a step voltage across the piezoelectric actuator and monitor- 
ing the velocity decay at a point on the beam structure. A Polytec Laser 
Scanning Vibrometer (PSV-300) was used to monitor the velocity decay at 
a point located on the surface of the beam. A dSPACE DS1103 system was 
used to generate the voltage step and capture the time response data from the 
vibrometer. The sampling frequency was set at 20 kHz. Within the dSPACE, 
a low-pass digital filter was used to remove high order structural modes and 
noise. The cut-off frequency for the digital filter was set at 230 Hz. After sam- 
pling the time response, the data was filtered with a high-pass digital filter 
to remove low frequency noise. The cut-off frequency was chosen to be 30 Hz 
so as to remove the dynamics of the first structural mode, which was not 
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Figure 4.16. Magnitude frequency response of the beam measured from an applied 
disturbance voltage to the resulting displacement. ( — ) Open- loop, - - shunted 





time (s) 

Figure 4.17. Transient step response of the beam 



targeted for damping. Results from the transient time regime are shown in 
Figure 4.17. 
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Feedback Structure of Piezoelectric Shunt 
Damping Systems 



5.1 Introduction 

The previous chapter introduced the concept of piezoelectric shunt damping 
for vibration control of flexible structures. From the discussions, it should be 
clear that most of the methods proposed in the literature for the design of 
shunt circuits are based on ad hoc procedures. It turns out, however, that 
the problem of piezoelectric shunt damping can be interpreted as a feedback 
control problem. The feedback structure associated with shunted piezoelectric 
transducers was first reported in [140, 139]. 

Having made the observation that there is an underlying feedback struc- 
ture associated with piezoelectric shunt damping we can propose new shunt 
impedances that are capable of damping structural vibrations in an efficient 
way. Furthermore, this can be done in a systematic way using the identified 
equivalence between the problems of piezoelectric shunt damping and collo- 
cated feedback control with an identical piezoelectric actuator and sensor. 

Two classes of resonant shunts are introduced in this chapter. These shunts 
are closely related to the two resonant controller structures introduced in 
Chapter 3. In particular, we will illustrate that the closed-loop system consist- 
ing of a piezoelectric transducer shunted by a resonant impedance is equivalent 
to the closed-loop system associated with an identical structure consisting of a 
piezoelectric actuator/sensor pair and a resonant controller. A direct implica- 
tion of this observation is that, in certain cases, it may be possible to replace 
a feedback controlled system by a shunted system, and thereby, reduce the 
number of necessary transducers by half. 

In addition to a thorough analysis of shunt damping systems, we will 
demonstrate the effectiveness of the proposed shunt impedances by imple- 
menting them on two piezoelectric laminated structures. 
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Figure 5.1. (a) A piezoelectric laminate shunted to an impedance Z(s); (b) Elec- 
trical equivalent of the same system 



5.2 Feedback Interpretations 

Consider the shunted piezoelectric transducer in Figure 5.1 (a) and its elec- 
trical equivalent in Figure 5.1 (b). It should be clear that compared to active 
control methods shunting the piezoelectric transducer with an impedance Z 
removes the need for an additional sensor. It will be shown, however, that this 
is achieved at the expense of having to deal with a slightly more complicated 
feedback control problem. 

To visualize the underlying feedback control structure, we need to identify 
a number of variables such as the control signal, the measurement, the dis- 
turbance and the physical variable that is to be regulated. Furthermore, we 
have to choose either Z(s) or Y(s) as the controller. The feedback structure 
can be identified by noticing that the current flowing out of the piezoelectric 
transducer may be written as: 

i(s) = - (v p (s) +v(s))C p s. (5.1) 

Furthermore, 

v(s) = Z(s)i(s). (5.2) 

Also, given the linear behavior of the system, we may write: 

Gyy^S^V^S^ T Gy W (s)w(s). (5.3) 

Here, G vw represents the transfer function from the disturbance to the 
induced piezoelectric voltage v p , when the piezoelectric transducer is short 
circuited. Also, G vv is the collocated transfer function of the piezo. That is, 
the transfer function from a voltage applied across the conducting electrodes 
of the transducer to the induced piezoelectric voltage v p . Furthermore, as 
explained in Chapter 2, typically 
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(a) 




(b) 



Figure 5.2. The feedback structure associated with the shunt damping problem in 
Figure 5.1: (a) Z(s) functions as the controller; (b) T(s) serves as the controller 



M 



G vv (s) = J2 






7 i 

s 2 + 2 QuJiS + uj, 



2 ’ 



and G vw has a similar structure, i.e. 

M 

c ( s) = Y ' ^ 

j f^s^ + ZCiUJiS + ujl 



(5.4) 



(5.5) 



G vv represents a collocated transfer function, hence 7* > 0 for i = 1, 2, . . . , M. 
However, due to the non-collocated nature of G vw , a similar assertion cannot 
be made, i.e. in general r]i may be non-positive. 

Equations (5.3), (5.1) and (5.2) amount to the feedback structure depicted 
in Figure 5.2(a). The block diagram suggests a slightly complicated feedback 
structure as the controller, Z(s) is itself inside an inner feedback loop. 
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Figure 5.3. Equivalent feedback loop 



If Y ( s ) is to be viewed as the controller, the block diagram of Figure 5.2(a) 
should be redrawn as in Figure 5.2(b). There are specific reasons as to why 
Z(s) or Y(s) should be chosen as a controller. Some of these reasons will be 
clarified in this chapter. 

The reader may notice that the feedback systems depicted in Figure 5.2 
(a) and (b) are indeed equivalent to the feedback control problem associated 
with the collocated system G vv (s), and a feedback controller K(s) depicted 
in Figure 5.3, where 



or alternatively 



sCpZ(s) 

1 T sCpZ(s) 


(5.6) 


1 

“ i + ivr 

1 + sC„ 


(5.7) 



An implication of this observation is that the shunted system in Fig- 
ure 5.4(b) and the feedback system in Figure 5.4(a) are essentially equivalent 
as long as K(s) and Z(s) are related according to (5.6) or (5.7). This observa- 
tion, however, could be misleading as it may lead the reader to the conclusion 
that having designed a controller for the former system, one may obtain an 
impedance for the latter. While this may be true in certain cases, such a pro- 
cedure may result in an impedance, or an admittance transfer function that 
is not implementable digitally. Therefore, more practical impedance design 
methods are needed. A number of techniques are discussed in the following 
sections. 

Identification of the underlying feedback structure associated with shunt 
damping is an important step in designing high-performance impedance 
shunts. In particular, the knowledge of this feedback structure enables us to 
address issues that would be very difficult to tackle otherwise. This includes 
problems such as: fundamental performance limitations in vibration damp- 
ing; dealing with actuator saturation; multivariable shunt design; robustness 
issues, etc. Some of these issues will be discussed in the remainder of this 
book. 
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K(s) 









(a) 




Figure 5.4. (a) Feedback control of a collocated system - notice the high-impedance 
buffer between the sensing patch and the controller, (b) A shunt damped system 
with a redundant open-circuited piezoelectric transducer. 



To this end, it should be mentioned that if the disturbance signal is applied 
to an identical collocated piezoelectric transducer, such as shown in Figure 5.5, 
then G vw = G vv . This configuration is frequently used in laboratory experi- 
ments. In fact, some experimental results obtained from a similar setup will 
be demonstrated at the end of this chapter. 

A point that needs to be clarified before this section is concluded is the 
effect of the redundant piezoelectric transducer in Figure 5.4(b). Existence of 
this patch is essential, if the systems depicted in two parts of this figure are 
to be equivalent. If the redundant transducer is removed, or short circuited, 
the transfer function G vv (s ) associated with Figure 5.4(b) should be replaced 
with G vv (s ), where 



G vv (s) 



G vv (s ) 

1 + G vv (s) 



(5.8) 



This should be of no particular concern since G vv (s ) represents a collocated 
system. Hence, any methodology developed for one system can still be ap- 
plied to the other. In most practical applications, of course, this redundant 
transducer should be removed. Indeed, a significant advantage of piezoelectric 
shunt damping, compared with using a collocated feedback control system, is 
that the need for a sensor is abolished. 



5.3 Feedback Structure of Passive Shunts 

In this section we attempt to apply the observation made in the previous 
section, on the feedback nature of piezoelectric shunt damping, to develop 
a deeper understanding of how efficient, or otherwise, certain passive shunts 
may be. 
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Figure 5.5. A structure with collocated piezoelectric transducers. 



The first shunt to be studied is a purely resistive shunt, that is, Z(s) = R. 
This amounts to the feedback system consisting of G vv (s ) and the feedback 
controller 



K r 



s 

s T (jJ c 



(5.9) 



where 



uj c 



1 

~RC~ P 



It can be observed that in this case the equivalent feedback controller is 
a first order high-pass filter. If the cut-off frequency of this filter, cj c , is cho- 
sen to be high enough to include several modes of G vv , the controller will 
approximate a differentiator within that bandwidth. In Section 3.4 it was 
demonstrated that a collocated system can be effectively damped using a dif- 
ferentiator. However, given the fixed nature of Kr , the gain that is needed to 
provide sufficient damping cannot be supplied. Consequently, a purely resis- 
tive shunt cannot act as an effective damper. Another possibility, of course, 
is to lower the cut-off frequency, by choosing a large enough resistor. In this 
case, the controller will, effectively, perform as a gain, which cannot provide 
sufficient damping either. 

If the shunt is a purely capacitive load, i.e. Z(s) = ^ , the equivalent 
feedback controller will be 



K(s) 



Cp 

C p + C‘ 



Therefore, the controller will be a fixed gain. This gain will always be less 
than one. Consequently, a significant level of damping cannot be expected 
from this shunt. 

A purely inductive load, Z(s) = Ls , results in the feedback controller 



K l 



S 2 + ’ 
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where 

1 

UJ '£ — , 

This controller has both its poles on the joo axis and is very ineffective in 
adding damping to the structure. Suppose the resonant system is made up of 
second order systems with no damping. That is, all poles of G vv are on the 
j uj axis. Then, all the poles of the closed- loop system will be on the juo axis 
too. Therefore, Kl adds no damping to the system. In reality, however, each 
mode of G vv has a small amount of damping. Introduction of this controller 
“may” add very small damping to the system. The added damping, however, 
will not be significant. 

Now, if the shunt consists of the series connection of a resistor i?, and an 
inductor L, the equivalent feedback controller will be: 



s(s + 2d t u e ) 

RL s 2 + 2 dgujgs + uj 2 ’ 



(5.10) 



where 






1 

777 



and 



de = 



R 

J 




It can be observed that the equivalent feedback control problem consists of 
G vv (s) and the controller K s RL (s ) which is a narrow band-pass filter centered 
around the resonance frequency of G vv (s). Hence, the controller would 
apply a high gain at that specific frequency, which would result in addition of 
extra damping to the structure, associated with that specific mode. The key 
difference between Krl and Kl is dg. This parameter can be manipulated in 
a way that the damping added to the structure is maximized 1 . 

For a shunt which consists of the parallel connection of a resistor R with 
an inductor L, the equivalent feedback controller will be 



Kr L ( s ) ~ 2 



s 2 T 2 d£(jj£S T (dg 



where uog is as above, and 




Given the structure of the equivalent controller, it is natural to expect this 
shunt to be equally effective in adding damping to the structure. 



1 This amounts to adjusting the resistor R 
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5.4 Reduction of Inductance Requirements by Adding a 
Parallel Capacitor 

It was noted in Section 4.2.3 that one way of reducing the size of inductors 
needed in piezoelectric shunts is to add a capacitor across the shunted piezo- 
electric transducer. Although we have demonstrated that almost any passive 
circuit can be digitally implemented using either a synthetic impedance or a 
synthetic admittance arrangement, in certain applications it may be benefi- 
cial to reduce the size of the necessary inductors to the point that they are 
commercially available. The discussion in this section targets this particular 
set of applications. 

Understanding the feedback structure associated with shunt damping sys- 
tems enables us to explain how a capacitor should be chosen, and what com- 
promises need to be made in order to make the system work properly. 

Consider the shunted piezoelectric transducer in Figure 5.6 (a), where a 
capacitor, C a dd , has been added in parallel to the shunt Z(s). It can be proved 
that the shunted transfer function from w to v v is [67]: 



L VpW 



(s) = 



C VW (<s) 






(5.11) 



where 



K(s) 

Ct 



aC T Z(s) 

i + c T z( s y 

Cadd T Cp, 



and 



a = 



Cp_ 

Ct 



Therefore, the system illustrated in Figure 5.6 (a) is equivalent to that in 
Figure 5.6 (b). According to the choice of Ct, the structure of the system is 
changed. Hence, this may adversely affect the performance obtainable from 
the system. However, it could have the advantageous property of reducing the 
size of necessary inductors. The following analysis of an RL shunt will clarify 
this point. 

If the piezo is shunted by a series RL circuit, i.e. if Z(s) = R + Ls, then 
the equivalent feedback controller is K S RL , given in (5.10). Adding a capacitor 
C a dd in parallel to this impedance, assuming that the resistive and inductive 
components are changed to R and L respectively, results in the equivalent 
controller 



~ s s(s + 2d t U> t ) 

RL s 2 + 2 d^ujis + o;| ’ 



(5.12) 



where 



1 
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Figure 5.6. (a) A shunted piezoelectric transducer with a parallel capacitor and 
(b) its equivalent circuit 



Table 5.1. Shunt circuit component values 



and 





Natural PZT 


Modified Cap. 


c P 


104.8 nF 


104.8 nF 


Cadci 


0 nF 


212.2 nF 


C 2 


104.8 nF 


317.0 nF 


L i 


41.8 H 


13.9 H 


l 2 


20.8 H 


6.9 H 


l 2 


41.8 H 


13.9 H 


R i 


1543 n 


514 n 


R2 


1145 R 


381 R 


di ~ 

2 V 


[Of 

L ' 



The original feedback controller can be retained by enforcing the following 
conditions: 



L = L — 

L C t 

= La 



and 



R 



R C t 

Ra. 



Thus, to ensure K RL (s ) = K s RL (s ), the inductive and resistive elements should 
be reduced by the same factor that the capacitance is increased. 

For multi-mode shunt circuits, such as those described in Chapter 4, similar 
results can be obtained. By increasing the effective piezoelectric capacitance, 
the resistance and inductance of each branch should be reduced. 
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Therefore, adding a capacitor in parallel to a piezoelectric transducer in- 
creases the capacitive part of the transducer. Inductive and resistive param- 
eters of the shunt can be manipulated in a way that the original feedback 
controller is recovered. This way, however, the original performance cannot 
be achieved since parameter a is always smaller than one. That is, the supple- 
mentary capacitance has the effect of reducing gain of the feedback controller. 

The above analysis is best demonstrated on a two- mode series shunt, sim- 
ilar to that described in Section 4.2.4. A summary of the component values 
is given in Table 5.1, where the first column represents the original shunt pa- 
rameters. By attaching a capacitance of C a dd = 212.2 nF across the shunt, 
a threefold reduction of component values can be achieved, as demonstrated 
in Column 2 of Table 5.1. Figure 5.7 shows the frequency response of the 
system for three cases: open- loop, shunted and shunted with a supplementary 
capacitance. Addition of this supplementary capacitance has resulted in a de- 
crease in the performance of the shunted system, as expected from the above 
analysis. 




Figure 5.7. Displacement magnitude frequency response of the shunted beam when 
the piezo is open circuited, shunted and a supplementary capacitance is placed across 
the shunt 







